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Abstract. We consider natural complex Hamiltonian systems with n degrees of freedom given 
by a Hamiltonian fimction which is a sum of the standard kinetic energy and a homogeneous 
polynomial potential V of degree k > 2. The well known Morales-Ramis theorem gives the strongest 
known necessary conditions for the Liouville integrability of such systems. It states that for each 
k there exists an explicitly known infinite set Mj- C Q such that if the system is integrable, then 
all eigenvalues of the Hessian matrix V"{d) calculated at a non-zero d E C" satisfying V'{d) = d, 
belong to Mj^. 

The aim of this paper is, among others, to sharpen this result. Under certain genericity assump- 
tion concerning V we prove the following fact. For each k and n there exists a finite set 3„ C Myt 
such that if the system is integrable, then all eigenvalues of the Hessian matrix V" [d) belong to 3,, 
We give an algorithm which allows to find sets 3„ jt- 

We applied this results for the case n = k = 3 and we found all integrable potentials satisfying 
the genericity assumption. Among them several are new and they are integrable in a highly non- 
trivial way. We found three potentials for which the additional first integrals are of degree 4 and 6 
with respect to the momenta. 
MSC2000 numbers: 37J30, 70H07, 37J35, 34M35. 
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1 Introduction 

The problem of the integrability of low-dimensional Hamiltonian systems has been a very 
active field of research in the last half century. Let us recall that a Hamiltonian system with 
n degrees of freedom defined by Hamilton's function H = H{q, p), where q = {qi, . . . , cj„) 
and p = (pi, . . . , p„), is integrable in the Liouville sense if the canonical equations 

d dH d dH . ^ 

admit n functionally independent and commuting first integrals. Investigations of physi- 
cally important systems and more or less systematic analysis of systems depending on a 
certain number of parameters gave rise to several new integrable systems. Among the tools 
used for such investigations two are most popular: the direct method and the Painleve test. 
In spite of their great successes, these techniques are rather limited and are not appropriate 



for a more ambitious program which can be formulated in the following way: 

Give a complete list of integrable Hamiltonian systems of a given class. In other words: formulate 

necessary and sufficient conditions for the integrability of a given class of Hamiltonian systems. 

Such a program could be attractive if the chosen class contains, on the one hand, phys- 
ically interesting systems, and on the other hand, it is large enough, and moreover is, in 
some sense, naturally distinguished. In this paper we consider natural Hamiltonian sys- 
tems given by the following Hamiltonian 

H=ltp' + ^i'j)' (1-2) 

^ ! = 1 

where V{q) is a homogeneous polynomial of degree k > 2. The canonical equations corre- 
sponding to the above Hamiltonian have the form 

5., = ^.^ = -V'(,), (1.3) 

where V'{q) denotes the gradient of V{q). We say that a potential V is integrable if the 
above canonical equations are integrable in the Liouville sense. 

The class of systems given by (|1.2|) satisfies the imposed requirements. Among homoge- 
neous potentials there are many important systems which come from physics, astronomy 
and other natural sciences, see examples in paper [2] and references therein. Further- 
more, thanks to works we know that if a non-homogeneous polynomial potential 
y = ^min + ■ ■ ■ + V^max/ where Vmin and Vmax are homogeneous terms of the lowest and the 
highest degree, respectively, is integrable, then potentials Vmin and ^max are also integrable. 
Thus the integrable homogeneous potentials can be considered as building blocks of more 
complicated non-homogeneous integrable models. 

Obviously, having in mind the described program, we must be equipped with tools 
strong enough to guarantee its realisation. Fortunately, we have at our disposal two meth- 
ods: the Ziglin [39; 40] and Morales-Ramis [27; 23] theories which are very strong and 
effective. In both theories obstructions to the integrability are derived from an analysis 
of the variational equations along a non-equilibrium particular solution of the considered 
system. In this paper we will use the Morales-Ramis theory that can be considered as alge- 
braic extension of the Ziglin method and its detailed description with many examples one 
can find in books see also iSSIIiliij. 



The main theorem of the Morales-Ramis theory is following. 

Theorem 1.1 (Morales,1999). If a complex Hamiltonian system is integrable in the Liouville sense 
in a neighbourhood of a particular non-equilibrium solution, then the identity component of the 
differential Galois group of the variational equations along this solution is Abelian. 

For a proof, examples and the discussion, see |l23i] . 

Remark 1.1 Jn the above theorem we have to distinguish two cases. If the variational equations 
are Fuchsian, or its irregular singularities correspond to finite points in the phase space, then we 
assume the integrability with first integrals which are meromorphic in a connected neighbourhood 
of the phase curve. If the variational equations have an irregidar singidarity corresponding to a 
point in the phase space lying at the infinity, then we assume the integrability with first integrals 
which are meromorphic in a connected neighbourhood of the phase curve and additionally have a 
meromorphic growth at the infinity. In the case of Hamiltonian systems considered in this paper it 
means that first integrals are rational. 
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Here we underline that the above theorem is very well suited to a study of the integra- 
bility of homogeneous potentials, because system il.3} admits, at least in a generic case, 
a certain number of particular solutions. For general Hamiltonian systems il.l} it is not 
completely obvious how to find particular solutions in order to have a chance to apply the 
Ziglin or Morales-Ramis method. 

1271] Morales and Ramis obtained very nice general integrability obstructions for 



homogeneous potentials with an arbitrary number n > 2 of degrees of freedom and an 
arbitrary degree of homogeneity k. As it was observed by H. Yoshida, system (|1.3|l has a 
particular solution of the form 

q{t) = cp{t)d, p{t) = q>{t)d, (1.4) 

provided that the scalar function (p[t) satisfies the following equation 

(p = 

and d G C" is a solution of polynomial equations 

V'{d)=d. (1.5) 

A vector d satisfying the above equations is called the proper Darboux point of the poten- 
tial. 

The variational equations along solution \1A} have the form 

x = y, y = -cp{tf-^V"{d)x. (1.6) 

Hessian ^"(rf) of the potential V calculated at a proper Darboux point d is a symmetric 
matrix. Let us assume that V'\d) is diagonalisable. Then there exists a complex orthogonal 
n X n matrix A such that the canonical transformation 

x = At], y = A^, 

transforms system (|1.6|) to the form 

Tii = ^h = -Ai(p{tf-^f]i, i = l,...,n, 

or simply 

iii = -\,cp{tf-\, i = l,...,n, (1.7) 

where (Ai, . . . , A„) are eigenvalues of V'\d). It can be shown that if the identity component 
of the differential Galois group of the above system is Abelian, then the identity component 
of the differential Galois group of each equation in the system is Abelian. 

Remark 1.2 As it was explained in the assumption that V'\d) is diagonalisable is irrelevant. 
That is, the necessary conditions for the integrability are the same: if the potential is integrable, 
then for each A G spectr F"(ti) the identity component of the differential Galois group of equation 
ij = —A(p{t)^^^f], is Abelian. 

For a given energy e G C* the phase curve Fg associated with the particular solution 
(p{t) is, for A: = 3,4 an elliptic, or for fc > 4 a hyperelliptic curve given by 

<p2 = ^ - /) , , = le. (1.8) 
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As it was observed by H. Yoshida |37l]. each of equations (|1.7|) can be transformed to the 
hypergeometric equation. It can be done by the following change of the independent 
variable 

t — >z:=^(p{tf. (1.9) 



After this transformation equations (|1.7|) read 



z(i-z)//;' + 



k-l 3k -2 



2k 



(1.10) 



where i = 1, . . . ,n. Equation (|1.10|l is a special case of the hypergeometric differential equa- 
tion, for which, thanks to works of Schwarz js^, Kimura fisl l and others, the differential 
Galois group is well known. 

This fact combined with Theorem 11.11 has allowed J. J. Morales-Ruiz and J. P. Ramis 
to formulate in |27l ] a general theorem concerning the integrability of Hamiltonian sys- 
tems with a homogeneous potential. Here, we formulate this theorem for a polynomial 
homogeneous potential. 

Theorem 1.2. If Hamiltonian system (|1.3|) with polynomial homogeneous potential V{q) of degree 
k > 2 is meromorphically integrable in the Liouville sense, then for a proper Darboux point the 
values of {k, Aj) for i = 1, . . . ,n belong to the following list 



1. 
3. 
5. 
7. 
9. 



5,- 



40 10 ^ ^' 



where p is an integer. 



k-l 



+ p{p + l)k 



(1.11) 



Remark 1.3 By Remark [I72l Theorem [172] is true independently of the fact whether matrix V"{d) 
is diagonalisable or not. However, ifV"{d) is not diagonalisable, new obstacles for the integrability 
appear. Namely, if the Jordan form ofV"{d) has a block 



/3(A) := 



A 1 0' 
A 1 
A 



then the system is not integrable. Moreover, if the Jordan form ofV"{d) has a two dimensional block 
/2(A), and A belongs to the first item of table (|1.11|) , then the system is not integrable. This fact was 
proved in jfl^/ . 

We denote by Mjt a subset of rational numbers A specified by the table in the above 
theorem for a given k, e.g., for fc > 5 we have 



P + -p{p-i) I pez [>u 



+ p{p + l)k 



peZ 



(1.12) 
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Let us make some remarks about practical aspects of applications of Theorem 11.21 for a 
study of the integrability of a potential depending on several parameters. At first, we have 
to find a proper Darboux point, i.e., to solve nonlinear equations l|1.5|). And just at this 
very initial step we meet a very fundamental problem. Namely, even if equations (|1.5|) 
have solutions, in most cases we cannot find their explicit form. Unfortunately, it is not 
the only problem that we have. Even if we find the explicit form of a Darboux point d, 
its coordinates depend on parameters in a very complicated way and, moreover, we obtain 
much more intricate expressions for eigenvalues A, of V"[d). So, it is at least doubtful 
if we can extract useful information from the restriction A, S M/^, and this is not only a 
technical problem. In fact, a restriction of the form A, G Mjt gives rise to an infinite number 
of potentials which form co-dimension one families in the parameters space. For each of 
these families, the necessary conditions for the integrability are satisfied, so we need to 
apply stronger theorems to check if within these families there are integrable potentials. 
Hence, it seems that we have to restrict applications of Theorem [L2] to very limited classes 
of homogeneous potentials. 

Our aim in this paper is to show how to overcome the above mentioned difficulties. The 
basic observation is following. The existence of a proper Darboux point together with the 
necessary conditions given by Theorem [L2] put restrictions on the potential no matter if we 
know explicitly the proper Darboux point or not. The question is, if we can deduce from 
this fact a kind of a global restriction on V. We give a positive answer to this question. But 
the most amazing fact is that this answer implies the following statement which is valid 
for all homogeneous polynomial potentials of degree k > 2 which satisfy certain genericity 
conditions: 

For each n and k there exists a finite set J^j^ C such that if V is integrable, then for each 
Darboux point d of V all eigenvalues of V"{d) belong to J^j^. Moreover there is a constructive 
algorithm which allows to find sets J„^k- 

The described general idea has found successful applications for systems with two 
degrees of freedom, see lHH; S 112] ■ So nowadays all integrable homogeneous potentials 
of degree 3 and 4 are knowrQl In this paper we show that the above idea has a nice 
generalisation for systems with n > 2 degrees of freedom. Moreover, we apply the obtained 
results to a study of the integrability of homogeneous potentials with three degrees of 
freedom of the homogeneity degree equal to 3. 

The methods applied in this paper force us to work with complex Hamiltonian systems. 
So, we assume that our system given by Hamiltonian (|1.2|) is defined on a complex linear 
symplectic space {C^",to) equipped with the canonical symplectic form 

n 

CO = Y^dqiAdpi. (1.13) 

i=l 

Moreover, we assume that potential V is also complex, i.e., V G C)t[^], where Cj^lq] denotes 
C-linear space of homogeneous polynomials of degree k. 

In order to make a reasonable classification, we have to divide all considered potentials 



into equivalent classes. To this end we proceed as in 12211 . Namely, let PO(n, C) be the 
complex projective orthogonal subgroup of GL(n, C), i.e., 

PO(n,C) = {A e GL(n,C), | AA^ = aE„, cc G C*}, (1.14) 

where £„ is n-dimensional identity matrix. We say that V and V are equivalent if there 
exists A G PO(n, C) such that V{q) = Va{(]) := V{Aq). Later a potential means a class of 
equivalent potentials in the above sense. 



^For k = i there exists a family of potentials depending on one discrete parameter and it is unknown if this 
family contains an integrable case or not, see L221 for details. 
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2 Darboux points 



In the Morales-Ramis Theorem 11.21 the considered particular solution has a very special 
form (|1.4|l . To find it, we have to solve a system of polynomial equations V'{d) = d. In 
our considerations it is important to take into account all particular solutions. An analysis 
of this problem leads us to a definition of a geometrical notion of a Darboux point of a 
homogeneous potential. It was introduced in [22] for systems with two degrees of freedom. 
However, the passage from two to three and more degrees of freedom introduces essential 
difficulties which force us to modify slightly this concept. In this section we give basic 
definitions and prove several general facts concerning Darboux points. 

Let us remark that the definition of a Darboux point given in this section is closeh 



closely 

illii 



related to the notion of the proper direction (radial orbits) introduced by Guillot lU 
for polynomial differential equations with homogeneous right hand sides. 

2.1 Basic notions of algebraic geometry 

At first we recall very basic notions from the algebraic geometry which are necessary for 
our further considerations. The main purpose is just to fix the notation. We refer the reader 



to [34l | for a more general and formal exposition. 

The set of common zeros of polynomials /i, . . . ,/s € C[x] := C[xi, . . is denoted by 
V(/i, . . . ,/s). It is called the affine algebraic set. An algebraic set is a finite sum of disjoint 
connected components. We say that point p G X := V(/i, ... ,/s) is a simple point of X, iff 

-ip) 



rank 



9/1 9/1 



dxi dx„ 



ra-dimXp, (2.1) 



where Xp is the component of X containing p. In particular, if X is a finite set, then p ^ X 

is a simple point, iff the above Jacobi matrix has the maximal rank. 

A polynomial / G C[x] can be uniquely written as a sum of homogeneous terms. The 

highest degree homogeneous term of / is denoted by /+ . 

We say that /i,. . .,/m G C[x] do not intersect at the infinity, iff the only solution of 
= ■■■=/+ = is X = 0. The following fact is well known, see e.g. l|9l: [l9l:l3ll ]. 

Proposition 2.1. If fi, . . . ,fm G C[x] do not intersect at the infinity, then set V(/i,...,/m) is 
finite. 

A point in the m dimensional complex projective space CP'" is specified by its homoge- 
neous coordinates [z] = [zq : • • • : z^], where z = (zq, . . . ,Zm) G C"'^^. If 

Ui := {[zo : ■ ■ ■ : z„,] G CP'" | z; ^ } for i = 0,...,m, (2.2) 

then 

m 

CP'" = IJ Ui, (2.3) 

and we have natural coordinate maps 

0i : CP-^' D Ui ^ C", e,i[z]) = ixi,...,Xm), 

where 

(XI x,„) = f^ 'f-). (lA) 

\ Z; Z; Z; Z; / 
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Each Ui is homeomorphic to C". It is easy to check that charts {Ui,9i), i = 0,...,m form 
an atlas which makes CP*" an holomorphic m-dimensional manifold. It is customary to 
choose one Ui, e.g., Uq, and call it the affine part of CP"". Then we define the hyperplane 
at the infinity 

Hoc := {[zo : • • • : z^] e CP'" I zo = 0} . (2.5) 
It is clear that CP'" = Uq U Hoo and n Hoo = 0. The following map 

0CO : Hoo ^ CP'"-^ 0eo([O:zi : ••• :z„]) = [zi : ••• :z„], (2.6) 

is the bijection, so 0oo(Hoo) = CP'"~^. The above consideration shows that CP'" is homeo- 
morphic to C" U CP'«-i. 

Let Fi,...,Fs e C[2] be homogeneous polynomials. Then their common zero locus in 
CP'", i.e., 

V(Fi,...,F,) :={[z]G CP'" I F,-(z)=0, 2 = 1,. ..,s}, (2.7) 

is called the projective algebraic set. The affine part of X := V(fi, . . . ,fs) is, by definition 
X n L?o. It is homeomorphic to the algebraic set ^(/i, . . . ,/s), where fi is a dehomogenisa- 
tion of Fi, that is 

fi{xi,...,Xm) =Fi{l,xi,...,Xm) for i = l,...,s. (2.8) 

Let us show this. Notice that if f G C[z] is a homogeneous polynomial, then F{z) is not a 
function of point [z] G CP"". In order to describe the hypersurface V{F) G CP'" as a zero 
set of a certain function defined on a subset of CP"', we define on Uq the following function 

F:Uo^C, [zo : • • • : z«] ^ F{[z]) := z^'^'^^Fiz), (2.9) 

and then 

V{F)nUo = {[z] G Uo I F([z]) =0}. 
The function F is called the function defining V(F) on Uq- If we set 

/ = F o : C"' ^ C, f{xi, ...,Xn):= F(l, xi,..., Xm), 

then 

0o(V(F)nLZo)=V(/). 

Poljmomial / is called the dehomogenisation of F with respect to the first variable. The 
above shows that 

0o(V(Fi,...,F,)nLlo) = V(/i,...,/,) 

as we wanted to show. 

Let [z] G Ui for a certain < z < n. We say that [z] is a simple point of V(Fi, . . . , Fg) iff 
0,([z]) is a simple point of V(/i, . . . ,/s), where fj is the dehomogenisation of Fy with respect 
to the z-th variable, for ; = 1, . . . , s. 

For the later use we will need explicit relations between representations of ■V(Fi, . . . , Fg) 
on two different charts. Let 

01 (V(Fi,...,F,)nL[i) = V(/ii,...,/i,), 

where hi is the dehomogenisation of F, with respect to the second variable. A relation 

between polynomials hi and can be found with a help of the transition map = QqoQ:^^. 
Namely, if 

eo{[zQ: ■ ■ ■ -.Zm]) = {xi,...,Xm), for [zo : • • • : Zm] G LTo, (2.10) 
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and 

0i([zo : ■ ■ ■ : z„,]) = (1/1,. . .,ym), for [zq : • • ■ : z^] G Ui, (2.11) 

then 

{x^,...,x^) = ^p{y^,...,y^) =(^y-l,...yjl\. (2.12) 

Vyi yi yi y 

Then, one can check that 

Hy) ■■= yT'^'MHy)) = yT'^'fi ^) , (2.13) 

Vyi yi yi / 

for i = 1,. . . ,s. 

In our considerations we invoke several times the Bezout Theorem. Let us recall its 
formulation, as the same name have several different theorems. 

Theorem 2.1 (Bezout). The number of solutions N of a system ofm homogeneous equations 

Fj{z) =0, i = l,...,m, 
in m + 1 unknowns z := (zq, . . . ,Zm) is either infinite, or is equal 

m 

N = YldesFi, 



i=i 



provided that solutions considered as points in CP*" are counted with their multiplicities. 

A precise definition of the multiplicity notion used in this theorem is highly technical 
and needs involved language of algebraic geometry. Here we refer the reader to Chap- 
ter 4 in [9]. For the purpose of this paper it is sufficient to know that a simple point has 
multiplicity one. 

2.2 General properties of Darboux points 

Let y be a homogeneous polynomial potential of degree k > 2, i.e., V ^ Cj^lq]. A direction, 
i.e., a non-zero d S C", is called a Darboux point of V iff the gradient V'{d) of y at d is 
parallel to d. Hence, d is a Darboux point of V iff 

dAV'{d)=0, d^O, (2.14) 

or 

V'{d) = 7d, d^ 0, (2.15) 

for a certain 7 G C. Obviously, if d satisfies one of the above conditions, then d = ocd 
for any a G C* satisfies them. However, we do not want to distinguish between d and d. 
Hence we consider a Darboux point d = {di, . . . ,d„) ^ C" as a point [d] := [di : ■ ■ ■ : dn] in 
the projective space CP"^^. 

The set D{V) C CP"^^ of all Darboux points of a potential V is a projective algebraic 
set. In fact, D{V) is the zero locus in CP"^^ of homogeneous polynomials Kg G Cj-[<7] 
which are components of ^ A V'{q), i.e. 

dV dV 

, := qiT: qi^:—, where \ <i < i <n. (2.16) 

'' oqj "dqi 

We say that a Darboux point [d] G D{V) is a proper Darboux point of V, iff V'{d) 7^ 0. 
The set of all proper Darboux points of V is denoted by V*{V). If [d] G V{V)\V*{V), 
then [d] is called an improper Darboux point of potential V. 
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We say that [d] is an isotropic Darboux point, iff 

dl + ---+dl = 0. (2.17) 

The set of all isotropic Darboux points of potential V is denoted by 1)o{V). 
We consider also two additional subsets of ^^"^^ 



Z(V):^V(V), S(V):=V(||^ ^) . (2.18) 

and we can prove the following relations. 
Proposition 2.2. For a non-zero y G C)c [q] we have 

1. V*{V) = V{V)\§{V), 

2. §{V) C Z{V). 

The first of the above statements follows just from the definition of T>*{V), and the 
second from the Euler identity for the homogeneous polynomial V. 

To describe the set T){V) in the affine coordinates we put, according to our convention, 

Ui := CF"-^ \ { e CP"-i I qi=0}, and 

ei:Ui^C"-\ i:= =0i([^]), (2.19) 

where 

Xi = ^, for / = l,...,n-l. (2.20) 

Lemma 2.1. On the affine chart {Ui,9i) we have 

9i{V{V)nUi) =V{g^,...,gn-i), (2.21) 
where polynomials gi, . . ■ ,gn-i £ C[x] are given by 

z;(x) := y(l,xi,...,x„_i), go:=kv-Y2xij—, (2.22) 

and 

dv 

gi •= 9^ - ^iSO' for z = 1, . . . , n - 1. (2.23) 

Moreover, [d] ^ D{V) n Ui is an improper Darboux iff its affine coordinates a := 0i([rf]) satisfy 
go{a) = 0. 

Proof. We know that 

D(y) = V(Ri,2,...,R„-i,„), (2.24) 

where Rj^j, for 1 < i < j < n are given by (|2.16|) . Hence, as it was explained in the previous 
subsection, we have 

0i(13(y)nLJi) = V(ri,2,...,r„_i,„), (2.25) 

where 

ri^j{xi,...,Xn-i) := Ri,j{l,xi,...,x„^i), for 1 < i < j < n. 
Direct calculations give 



q\-'go{x), (2.26) 
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and 



^^"^ 3^ = '?'r^ (^'-1 + ^'-1^0 )' 2 = 2, . . . , n. (2.27) 

Using the above formulae we easily find that 

= gi' for i = l,...,n-l, 

and 

= Xigj - Xjgi, for 1 <i < j <n-l. 

It shows that 

VK2,...,r„_i,„) = (2.28) 

and thus 

0i(D(y)nLJi) = V(^i,...,g„_i). (2.29) 

To prove the last statement of the lemma let us assume that [d] is a Darboux point such 
that [d] e Ui with affine coordinates a = 6i{[d]). Then, gi{a) = 0, for f = 1, . . . , n — 1. If 
[d] is an improper Darboux point, then all partial derivatives of V vanish at d. From H2.26|l 
it follows that if 

^(d)=d^igo(5)=0, 

then necessarily go{a) = 0, because di ^ 0. On the other hand, if gi{a) = 0, for / = 
0, ...,n — 1, and 0^^{a) = [a], with a := (l,fli,. . .,fl„_i), then from (|2.26|) and (|2.27|) we 
have ^'(fl) = 0. Since [a] = [d], this shows that [d] is an improper Darboux point. □ 

The construction which we are going to describe now plays a very important role in 
our considerations. Roughly speaking, the idea is to associate with a Darboux point which 
is located in CP"^^, a finite set of points in CP". This procedure is a kind of blowup. 

Let us define the following n homogeneous polynomials of n + 1 variables q := {qo, qi, . . . , 

Fi := — - q^-^qu i = l,...,n, (2.30) 

oqi 

and an algebraic set T>{V) = V{Fi, . . . , F„) C CP". 

Assume that [d] G T)*{V). Then there exists 7 G C*, such that V'{d) = jd, so k - 2 
points [ ^-^^ : di : ■ ■ ■ : dn] E CP" belong to T>{V)- These points are well defined as they do 
not depend on a representative for \d] . If [d] is an improper Darboux point, then it defines 
just one point [Q : di : ■ ■ ■ : dn] E CP" which is a point of I'(V). 

Set T>{V) is not empty because it contains point [do] := [1 : : ■ ■ ■ : 0]. If [d] = [rfo : ^^i : 
• • • : dn] G T){V)\ {[do]}, then [d] = [di : ■ ■ ■ : dn] is a Darboux point of V. Moreover, if 
do 7^ 0, then [d] is a proper Darboux point. 

The natural projection 

TT : CP" \ {[do]} ^CP"-\ Tiilqo : qi : ■ ■ ■ : qn]) = [qi : ■ ■ ■ : qn], (2.31) 

maps 1){V) \ {[do]} onto T){V), and the intersection of the inverse image 7r^^([d]) of a 
Darboux point [d] E T>{V) with T){V) is a finite set. We define also 

7r:S(y)\{[do]}^2)(n (2.32) 

putting 7r([d]) := 7r([d]) for [d] G T){V) \ {[do]}. That is, n is the restriction of n to 
T){V)\ {[do]}. This construction is illustrated in the Figure IZTl 
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Example 2.1 For potential 



(2.33) 



! = 1 



S(V) = 0, so D*(V) = I'(V). Let us calculate the number of points 0/5(1^). Since Fj 
qifor i = 1, . . . ,n, it is easy to find that card D(y) = {k — 1)". Hence, cardl'(y) 



,lc-l Jc-2 



cardD*(F) := D{n,k) where 



D{n,k) 



{k-iy- 



(2.34) 



Next lemmas show how to use the above construction for proving several facts which 
are not so obvious. 

At first we distinguish generic potentials. We say that potential V is generic iff all its 
Darboux points are proper and simple. This definition is justified be the following lemma. 

Lemma 2.2. The set of generic potentials 9n,k C Ck[q] of degree k is a non-empty open set in C)t[^]. 
A generic V G Cjc[<7] has D{n,k) proper Darboux points. 



Proof. We proceed as in the proof of Proposition 4 in |12l], see also [11]. Except for point 
[do], every point of l3(y) is a Darboux point. Moreover, every proper Darboux point gives 
k — 2 points of 23(V). Hence, if hypersurfaces V(F;) C CP" intersect at a finite number 
of points in CP" and all these intersections are transversal, then by the projective Bezout 
Theorem 12.11 we have (A; — 1)" of such intersections which give D{n,k) Darboux points. 
Example 12.11 shows that for n > 2 and A: > 2 set 5n,k is not empty. Fact that Sn,k is open 
in Ck [q] considered as a finite dimensional C-vector space with the standard topology, was 
proved in Jsl], Lemma 2. □ 



Let us note that this lemma is a generalisation of Corollary 2.1 in 1221 1. when the case 
n = 2 was considered and then a generic potential has D{2,k) = k proper Darboux points. 




n-l 



CP" I 



Figure 2.1: Sets V{V) C CP" and D{V) C 



n-l 



Lemma 2.3. For an arbitrary V G Ck[q] the set of Darboux points T>{V) is not empty. 
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Proof. lfD{V) =0, then T>{V) = {[do]}. Then, by the Bezout theorem the multiplicity of 
[do] is greater than one. However, as it is easy to notice, [do] is a simple point of I'(y). A 
contradiction finishes the proof. □ 

On the affine chart {Uo,9o), the affine part of 'D{V) is an algebraic set . . . ,/„), 
where /, is the dehomogenisation of f„ i.e., 

dV 

fi{cii,...,qn) ■■= Fi{l,c]i,...,q„) = cji, for i = l,...,n. (2.35) 

OCjj 

Notice that the polynomials /i, . . . ,/n intersect at the infinity iff the gradient V' vanishes at 
a certain d ^ 0. But such d gives an improper Darboux point [d]. On the other hand, the 
points of intersections of I'(y) with the hypersurface at the infinity Hoo correspond exactly 
to the improper Darboux points, more precisely 

e^{T){v)nHoo) = v{v)\v*{v). 

The lemma below gives a necessary condition for the finiteness of set D{V). 

Lemma 2.4. If a potential V does not have the improper Darboux points, then it has a finite number 
of Darboux points. 

Proof. Under assumption of the lemma we have Ti{V) n Hoo = 0- So, all points of D{V) 
are located in the affine part of CP", and their affine coordinates are points of V(/i, . . . ,/„). 
Since, by the assumption, fi,...,f„ do not intersect at the infinity, by Proposition 12. 1[ set 
^{fi, ■ • ■ is finite. Hence, T){V) as well as V{V) are finite. □ 

Remark 2.1 Let [d] € 'D{V) be a Darboux point of potential V. Then for each A G PO(n, C), the 
potential Va equivalent to V has Darboux point [Ad]. This fact allows us to fix coordinates of one 
Darboux point. Namely, if [d] is not isotropic, then we can assume that d = (0, ... ,0, 1). If[d] is 
isotropic, then we can assume that d = (0, . . . , 0, i, 1). 

For a Darboux point [d] G "^{V) we can calculate eigenvalues Ai(d),. . .,A„{d) of the 
Hessian matrix V"{d). However, numbers A;(d) are not well defined, as they depend on a 
representative d of the Darboux point [d] . There are several possibilities to define properly 
the quantities related to the eigenvalues of V"{d) which do not depend on a choice of a rep- 
resentative of [d] . However, because of some historical reasons and the convention widely 
accepted in the literature, we choose the one which is a simple normalisation. Namely, if it 
is not otherwise stated, we always assume the following 

Assumption 2.1 If [d] is a proper Darboux point, then the chosen representative is such that it 
satisfies V'{d) = d. 

If [d] is an improper Darboux point, then the representative of [d] can be chosen arbi- 
trarily. 

Before we formulate our next lemma let us introduce the following auxiliary system 

±-cj = f{q):=V'{q)-jq, (2.36) 

where 7 G C. If [d] is a Darboux point, then f{d) = for a certain 7, i.e., d is an 
equilibrium of (|2.36[) . If 7 7^ 0, then [d] is a proper Darboux point and for 7 = is an 
improper one. On the other hand, if d is an equilibrium of (|2.36|) for a certain 7, then [d] 
is a Darboux point. Under Assumption 12.11 if [d] is a proper Darboux point, then d is the 
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equilibrium of (|2.36|l with 7 = 1, and the Jacobi matrix of the right hand sides of (|2.36|) at 
d is f'{d) = V"{d) — En, so its eigenvalues are A,((i) := A,((i) — 1, for i = 1, . . . ,n. 

We make a change of coordinates in the system (|2.36|) . Namely, assuming that qi ^ 0, 
we put 

q^, Xi = ^ for / = l,...,n-l. (2.37) 



^0 



Then x := (xq, Xi,. . .,x„_i) are coordinates on C" \ G C" | qi = 0}. Notice that x :- 
(xi, . . .,x„_i) = 9i{[q]). In these coordinates system (|2.36|) reads 



^0 



woix) : 
'Wj{x) := x; 



7^0 





4 Vo. 



^i, for / = 1, 



1, 



(2.38) 



where go, . . ■ ,gn-i S C[x] are given by (|2.22|) and (|2.23[) . It is clear that the Jacobi matrix 
f'{x) has the same spectral properties as 



J{x) 



d{wo,. . .,W„-i) 

9(xo,...,x„_i) 



(x). 



(2.39) 



However, it appears that using }{x) we obtain simpler proofs of some facts. 

Under Assumption 12.11 we have the following simple but important consequences. 

Lemma 2.5. If [d] is a proper Darboux point of homogeneous potential V of degree k > 2, then 

q{t) := (p{t)d, p{t) := (p{t)d, (2.40) 

is a solution of Hamilton's equation (|1.3|) provided (p = —(p^^^. Moreover, V"{d) ■ d = A„d with 
An = k — 1, and if additionally [d] is isotropic, then A„ is a multiple eigenvalue ofV"{d). 
If [d] is an improper Darboux point, then 

q{t) := td, p{t) := d, (2.41) 

is a solution of Hamilton's equations \\3\ . Moreover, V"{d) ■ d = And, with An = 0, and if 
additionally [d] is isotropic, then A„ is a multiple eigenvalue ofV"{d). 

Proof. The fact that (|2.40|) and \2A1\ are solutions of (|1.3|) can be checked directly. Moreover, 
the fact that d is an eigenvector of V'\d) with the prescribed eigenvalue follows directly 
from the Euler identity. Hence, we have only to show that if [d\ is isotropic, then A„ is a 
multiple eigenvalue of y"(d). 

Let us assume that \d] is a proper and isotropic Darboux point. We show that matrix 
V"{d) — En has the multiple eigenvalue A„ = k — 2. Here we notice that this matrix is, 
under Assumption 12.11 exactly the Jacobi matrix of the right hand sides of the auxiliary 
system \2.36) calculated at equilibrium point d. Hence eigenvalues of V'\d) — E„ and their 
multiplicities are the same as eigenvalues of the Jacobi matrix of system H2.38|) calculated 
at an appropriate point. Direct calculations give 



-^ + {k-i)xl-^^ 

ik-2)xl-'gi 



L (fc-2)x^-3^„_i 



Jc-2 
"^0 



dlgl 



^0 



dn-lgl 



,1-2 



,1-2 



dn-lgn- 



(x). 



(2.42) 



At first we assume that [d] := [1 : i : : • • • : 0] is an isotropic Darboux point, so Xq = 7 = 1 
and d = 9-i{[d]) = (i, 0, . . . , 0) . Thus we have 



~ ~ dv ~ ^ ^ 

go{d) = kv{d) - ig^(t^) = 1/ giW = -ikv{d), 



(2.43) 
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and 



dv 



Si = ^(^)' for / = 2,...,n-l. 



(2.44) 



Since at a proper Darboux point we have gi{d) = for i = 1, . . . ,n — 1, we obtain also the 
following conditions 



v{d) = 0, 



dv 



(d) = 0, for / = 2, . . . , n - 1, 



as well as 



and finally 



go{d) 



Jid) 



'k-2 digo 
3x^1 



dig„_i 
But from conditions (|2.45|) it follows that 



dn-lgO 
d„-igi 



Jn-lgn-1 



id). 



(2.45) 



(2.46) 



(2.47) 



|l{3) = -(t-2)i||(3)=*-2. 



and 



(d)=0, for i = l,...,n-l. 



(2.48) 



Thus, in the second row of ]{d) all but the second element are zero, and the second element 
is A; — 2, so /c — 2 is an eigenvalue of /(d) with the multiplicity at least two. 

If [d] is an arbitrary proper and isotropic Darboux point, then we can always change 
the reference frame with the help of an element of PO(n, C) in such a way that in the new 
frame d = (l,i,0, . . . ,0) G 0"+^. 

If \d\ is an improper and isotropic Darboux point, then we can assume also that it is 
[d] = [1 : i : : • • • : 0], and then we proceed as in the case of the proper Darboux 
point. For an improper Darboux point 7 = 0, and gi{d) = for i = 0, ... ,n — 1, where 
d = ei{[d\) = (i,0,...,0). Thus we have 



~ 9y ~ 

v{d) = 0, and — {d) =0, for i = 1, . . . , n — 1, 



(2.49) 



and 



m 



digo 
9igi 



9)!-igo 

^n-lgl 



{d), 



(2.50) 



.0 ^Ign-l ■ ■ ■ d„_ig„_i_ 

as Xq = 1. But using (|2.49|) and definition of gi, one can easily find out that 

^(d)=0, for f = l,...,n-l, 
dXi 



so indeed is a multiple eigenvalue of J{d). 



□ 
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We underline that the presence of more than one eigenvalue equal to k — 2 does not 
imply automatically that the Darboux point is isotropic. 

Next we explain relations between the properties of a Darboux point [d] G T){V), and 
properties of points 7r^^([d]) C T){V). 

Lemma 2.6. Assume that set D{V) is finite. Then point [d] G T>*{V) is a simple point ofViV) 
iff each point [d] G 7r^^([d]) is a simple point ofT){V). 

Proof Let [d] G T)*{V) be a proper Darboux point satisfying V {d) = d. Without loss of 
generality we can assume that di ^ 0. Then the affine coordinates of [d] on the chart 
{Ui,di) are a := (fli, . . . ,fl„_i) = 9i{[d]) G C""\ where 



for f = 1, . . . , n — 1, 



(2.51) 



see formulae (|2.19|l and H2.20|) . Thus, \d\ is simple iff 

detg-(a).= det ^te-;;;-f;;-;> W/o. 



Now, let us take point \d\ := [1 : di 



Of course, [2] G n'^{\d\). 



Moreover, the affine coordinates of \d\ on the chart in CP" are d. On this chart 

%{y) is given by V(/i, ...,/„), see (12351) . Thus [3] is simple iff 



det/(d) = det(y"(<i) - £„) ^ 0. 



(2.52) 



Notice, that /(<7) coincides with the right hand sides of system l|2.36|l with 7 = 1. Thus, 
matrix f{(]) is similar to the matrix ]{x) defined by l|2.39|) . with x = (xq, . . . , given 
by (|2.37|) . For q = d, we have x = a := {ao, fli, . . . , fl^-i), where aq = di and fl, for i > are 
given by (|2.5Hl . Matrix J{x) is determined by (|2.42|l (with 7 = 1), so we have 



/(«) 



«0 ^3j!-1^0 



/c-2; 



^0 



(2.53) 



where gi = gi{a), for / = 0, . . . , n — 1), see !|2.38|l and below. But because [d] is a Darboux 
point, a is an equilibrium of system (|2.38|l , so gi{a) = ■ ■ ■ = gn_i{a) = 0, and «o^go(«) = 
1. Thus we obtain 



/(a) 



and 



(fc-2) al-%go 
det/'(d) = det/(a) = - 



«0 ^^n-lgO 

al'%i-igi 



k—2z\ 
% On-lgn-1. 



(2.54) 



(2.55) 



Thus we show that [d] is simple if and only if its one preimage, namely [d], is simple. Let 
[dj] := y : di : ■ ■ ■ : dn], for j = 0, . . . ,k — 3, 
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where e is a primitive {k — 2) -root of the unity. Then 

n-\[d]) = {[do],...,[dk-3]}. 
Because the affine coordinates of [dj] are dj := e^'d, we have 

f{dj) = V"{e-id) - En = £-^My"(d) - En = fid). 

Hence our conclusion is vaHd for an arbitrary element of the preimage of [d], and this 
finishes the proof. □ 

If [d] is an improper Darboux point, then the situation is more complicated. We can 
always assume that di 7^ 0. Then the affine coordinates of [d] are a = 9i{[d]), where a is 
given by (|2.5H) . Since [d] is improper Darboux point gi{a) = for i = 0, . . . ,n — 1, and 
thus, as direct calculations show, we have 

g'(a) = v"{a) ■ [En-i + a • fl^j , 



(2.56) 

where we consider a as a one colunm matrix. One can show that matrix a • is diag- 
onalisable, its one eigenvalue is A = + ■ • • + a^^, and the remaining ones are zero. 
Thus 

detg'Cd) = (1 + fl2 + . . . + fl2_^) detv"{a). (2.57) 

Note that 1 + fl^ + • • ■ + fl^^ = iff [d] is isotropic. Hence, we have the following. 

Proposition 2.3. Assume that set T>{V) is finite, and let [d] G §(^) be an improper Darboux 
point. Then it is a simple point ofT)(y) iff it is not isotropic and detc"(a) ^ 0. 

The preimage 7r^^([d]) of an improper Darboux point [d] = [di : ■ ■ ■ : d„] is just one 
point, namely [d] := [0 : d-[ : ■ ■ ■ : d„]. Again we assume that di ^ 0. Now, we cannot use 
chart (LZq/^o) in order to check whether [d] is simple. We will use chart {Ui,9i). Hence, on 
the chart (LZi, set D{V) is V{hi, . . ., hn) where polynomials h = [hi, . . . , hn), according 
to (I2l3l) . are 

\yi yi yiJ 

for / = !,...,«. Direct calculations show that 



(2.58) 



h,iy) := kv{y) - Ey,^(y) - y\-\ h{y) := ^(y) - y^y\-\ (2.59) 
for 2 < i < n, where y := (1/2/ • • - fyn), and 

v{y) := y(i,i/2,...,yn). 

With the above formulae, we can calculate the Jacobi matrix h'{y). It has the following 
form 



h'iy) 



-{k-2)y\ 



k-3 



fk-2)yry\-^ ''^{y)-y\-^ 



dyl 



dhi 
dy2 



(y) 



dp 



(y) 



{k-2)yny\ 



k-3 



dyndy: 



iy) 



dy2dy3 



dyndy 2 



iy) 
iy) 



dhi 
dyn 



iy) 



dyidyn 



iy) 



dyl 



iy) - v\ 



.k-2 



(2.60) 
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where 



for 2 < j < n. 



(2.61) 



Notice that = a := {0,a2, ■ ■ ■ ,a„), where 

di 

Ui = —, for i = 2,. . . , n. 
Moreover, because [d] is an improper Darboux point, we have 



-(a)=0, for ; 
where we denoted a := {aj, ... ,an). This gives 



= 2, 



,,n. 



Let A: = 3. Using the above formulae we find that 



h' {a) = 



-1 —v"{a) ■ a 



-02 



(5) 



v"{a) ■ a 
{a) 



—a 



"'d^'v' 



■al)detv"(a). 



(2.62) 



Hence, as it is easy to see, we have 

det/i'(fl) = - {l + al + ---+al)detv"{a). (2.63) 
The above calculations show the following. 

Proposition 2.4. Assume that k = 3. Then the preimage d := fc~^{[d]) of an improper Darboux 
point [d] e 2)(V) is a simple point ofT){V) iff [d] is a simple point ofD{V). 

For > 3, the point d is not a simple point of D(V), as all elements of the first column 
of matrix /i'(fl) vanish. For such a case we have the following. 

Proposition 2.5. Assume that k> 3. Then 

rankfi'(fl) = rank (a) < n — 1. 

Proof. For A: > 3 we have 



h'{a) 



-v"{a)-a 



z;"(a) • d 
(fl) 



dyidyn 







(«) 



(2.64) 



The first row of this matrix is a linear combination of the remaining ones, so our claim 
easily follows. □ 

With the above notation, we say that [d] G S(V) = 2)(V) \ T>*{V) is minimally degener- 
ated iff rank (a) = n — 1. 
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2.3 Obstruction for the integrability due to improper Darboux point 



As we have already explained, the proper Darboux points are important because they give 
particular solutions of the considered canonical equations (|1.3|) and, thanks to this, we have 
the necessary conditions for the integrability given by Theorem 11.21 The question is if we 
can obtain any obstruction for the integrability from the existence of an improper Darboux 
point. Lemma 1231 shows that the improper Darboux points give a particular solution (|2.41|l 
of the considered canonical equations (|1.3|) . However, this solution has an extremely simple 
form and one can doubt if using it we can obtain any obstruction for the integrability. In 
fact, it is easy to notice that the monodromy group of the variational equations along 
solution (|2.4Hl is trivial. Thus, in the frame of the Ziglin theory we do not obtain any 
obstacles for the integrability. However, situation is different if we are using the Morales- 
Ramis theory. The following theorem demonstrates the power of the Morales-Ramis theory, 
and it amazingly shows that the improper Darboux points give, in some sense, stronger 
obstacles for the integrability than the proper ones. 

Theorem 2.2. Assume that a homogeneous potential V G C^lq] of degree k > 2 admits an improper 
Darboux point [d] S CP"^^. IfV is integrable with rational first integrals, then matrix V"{d) is 
nilpotent, i.e., all its eigenvalues vanish. 

Proof. We prove the theorem by contradiction. We assume that there exists a non-zero A G 
spectrV"{d). Thus, in an appropriate base, the variational equations along solution (|2.41|) 
contain equation if = —Xt^^'^t]. As it was shown in [22] the differential Galois group of this 
equation is SL(2, C) and hence the identity component of the differential Galois group of 
all variational equations is not Abelian. Thus, by Morales-Ramis Theorem 1 1.1[ the system 
is not integrable. A contradiction finishes the proof. □ 



The above theorem gives a very strong generalisation of Theorem 2.4 in [221] for systems 
with an arbitrary number of degrees of freedom. 



3 Basic Theorems 
3.1 Formulation 

Let \d\ be a proper Darboux point of potential V . Then, thanks to Assumption 12.11 we have 
well defined eigenvalues ki{d),. . .,\„{d) of the Hessian matrix V"{d). According to our 
convention \„{d) = k — 1 is the trivial eigenvalue. Let \ {d) = (Ai(d), . . . , A„_i(d)). Hence 
we have the following mapping 

V*{V) 3 [d] I — > \{d) G C"-\ (3.1) 

Assume that 'D*{V) is finite. Then the image of D*(y) under the above map is a finite 
subset of C"^^. The question is if we can find a potential V of degree k such that the 
elements in the image have values prescribed in advance. We show that the answer to 
this question is negative. More precisely, we prove that among X{d) taken at all proper 
Darboux points [d] G D*(V') a certain number of universal relations exists. These relations 
play the fundamental and central role in our considerations. 

To formulate our first theorem we define A{d) = {Ai{d), . . . , A„_i(d)), where A,(d) := 
A;(d) — 1 for i = 1, . . . ,n — \. By r,- for < r < n — 1, we denote the elementary symmetric 
polynomials in (n — 1) variables of degree r, i.e., 

r 

Tr{x) ■.= Tr{xi,...,Xn-l) = ^ Yl^'s' 1 < r < M - 1, 

l<!i<---<;V<f!-ls=l 
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and To(x) := 1. 

Our first theorem gives the explicit form of the above mentioned relations among A{d), 
[d] ^ 1){V) for a generic potential V. 

Theorem 3.1. Let V S Ck[q] be a homogeneous potential of degree k > 2 and let all its Darboux 
points be proper and simple. Then 

,n-l. 

The proof of the above theorem is given in the next subsection. We prove at first 
the existence of relations i3.2^ and later we show that from them relations (|3.3|) follow 
immediately. In effect Theorem 13. 1 1 gives n independent relations. 

Let us explain the importance of Theorem 13. II To do this we need more definitions. 

Let S„, denote the set of all unordered tuples A = (Ai, . . . ,Am), where A,- G C for i = 
l,...,m. For M > 0, the symbol denotes the set of all unordered tuples (Ai, . . . , Am), 
where A, G Gm, for z = 1, . . . , M. 

We fix A: > 2 and n >2, and say that a tuple A S Qn-i is admissible iff A, = A, • + 1 € M]^ 
for / = 1, . . . , n — 1. In other words. A; is admissible iff A, + 1 belongs to items, appropriate 
for a given k, in the table of the Morales-Ramis Theorem ll.2[ for i = 1, . . . ,n — 1. We 
denote the set of all admissible tuples by A„^k. If the potential V is integrable, then for 
each [d] S D*{V), the tuple A{d) is admissible. The set of all admissible elements A„j( is 
countable but infinite. 

If the set of proper Darboux points of a potential V is non-empty, and N = card 'D*{V), 
then the N-tuples 

ii{V) := (A(d) I [d] e V*{V) ) G etv (3.4) 

is called the spectrum ofV. Let A^j^ be the subset of C^_i consisting of N-tuples (Ai, . . . , An), 
such that A,- is admissible, i.e.. A, S A^^i^, for i = 1, . . . ,N. We say that the spectrum ^{V) 
of a potential V is admissible iff XL (y) S A^j^. The Morales-Ramis Theorem 11.21 says that 
if potential V is integrable, then its spectrum ^{V) is admissible. However, the problem is 
that the set of admissible spectra A^i^ is infinite. We show that from Theorem 13. II it follows 
that, in fact, if V is integrable, then its spectrum i^{V) belongs to a certain finite subset J^j, 
of A^i^. We call this set distinguished one, and its elements distinguished spectra. 

Theorem 3.2. Let potential V satisfy assumptions of Theorem 13.11 If V is integrable, then there 
exists a finite subset 'i^^. C A^,^, where N = card 'D*{V), such that L{V) G 3^^. 

Informally speaking, for fixed k and n, we restrict the infinite number of possibilities in 
each line of the Morales-Ramis table to a finite set of choices. 

The importance of the above theorem immediately forces us to ask if it is possible 
to relax assumptions of Theorem 13.1! and preserves its conclusions. Below we give two 
generalisations of Theorem 13.11 In the first one we restrict ourselves to the case k = 3, but 
we do not assume that all Darboux points are proper. However, for this generalisation we 
obtain a smaller number of relations. 

Theorem 3.3. Let V S C3[<7] be a homogeneous potential of degree k = 3 and let all its Darboux 
points be simple. Then relations (|3.2|) and ^3.3} with r = 0, . . . ,n — 2, are satisfied. 



and 



for r = 0, . . . 
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In the second generalisation of Theorem 13.11 we do not assume that all Darboux points 
are proper and A: > 3 can be arbitrary. 

Theorem 3.4. Let V G C|c[^], be a homogeneous potential of degree k > 3 which satisfies the 
following conditions 

1. if[d] e T)*{V), then [d] is simple, 

2. if[d] S D(y) \ D*(y) then [d] is minimally degenerated. 

Then relations (|3.2|l and with r = 0, . . . ,n — 2, are satisfied. 

The fact that Theorems 13.31 and 13.41 give a smaller number of relations than Theorem 13. II 
does not weaken their strength. Namely, we show that the following result holds true. 

Theorem 3.5. Let potential V satisfy assumptions of either Theorem 13.31 or Theorem 13.41 If V 
is integrable, then there exists a finite subset jj^^ c A^j^, where N = cardD*(y), such that 



3.2 The Euler-Jacobi-Kronecker formula and multi-dimensional residues 

In proofs of Theorems 13.11 13.31 and 13.41 formulated in the previous subsection, we use the 
classical Euler-Jacobi-Kronecker formula and its known generalisations. In these generali- 
sations the multi-dimensional residues are used. We need to calculate these residues effec- 
tively. In this subsection we recall basic facts about the multi-dimensional residues and the 
Euler-Jacobi-Kronecker formula. For details the reader is refered to |fll:E4: 

Let fi : C" D U C, where U is an open neighbourhood of the origin, be holomorphic 
functions for i = l,...,n, and x = be an isolated common zero of /,. We consider 
differential n-form 

^ := . . , . dxi A • ■ ■ A dx„, (3.5) 

where p : 17 ^ C is a holomorphic function. The residue of the form a; at x = can be 
defined as 

res(a;,0) := ,^ ., / oo, (3.6) 

where 

Y:={xeU\ |/i(x)| =£i,...,|/„(x)| =£„}, (3.7) 
and El, ...,£,! are sufficiently small positive numbers. The orientation of T is fixed by 

d(arg/i) A--- Ad(arg/„) >0. (3.8) 

Let us denote / := (/i, . . . ,/n). It can be shown that if the Jacobian det/'(0) 7^ 0, then 

The following theorem gives the classical Euler-Jacobi-Kronecker formula, see e.g. HI]. 

Theorem 3.6 (Euler-Jacobi-Kronecker). Let fi,...,f„ S C[x] be non-constant polynomials such 
that V(/) := V(/i, . . . ,/„) is finite and all points of this set are simple. If f\, . . . ,f„ do not intersect 
at infinity, then for each p E C[z] such that 

degp<£deg/,-(n + l), (3.10) 

;=1 
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we have 

X: res(a;,d)= = 0" (^.ll) 

dGV(/) dGV(/) '^"Z 

The above theorem is not sufficient for our investigations. We have to consider cases 
when fi,...,fn have intersections at the infinity as well as cases when intersections of 
fi,...,fn are not simple. 

The homogenisations of /, are given by 

Fi{zo,zi,...,Zn) ■.= ZQ''^^'fi (—,..., — ] , for i = l,...,n. (3.12) 

They define the projective algebraic set "\7(F) := V(fi,. . C CP" whose affine part is 
homeomorphic to V(/). Next we extend the form a; to a rational form O defined on CP". 
To this end we consider co as the expression of Q on the chart ( Uq, Oq). In order to express 
O on other charts we use the standard coordinate transformations. For example, if we are 
using formulae from Section l2Tl on {Ui, 6i), the form Q is given by the puUback co := ^p*co. 
To write down co explicitly we put 

h,{y) := yT'^'MHy)) = vt^^'h ^ ^) . O-B) 

Vyi yi yi/ 



for i = \, . . . ,n, and 



Then we obtain 



where 



Av) ■.= yT''p{m)=yT''p (^.f ^) • (3.14) 

\yi yi yij 

^ = - I, ( ( ^ dyl ^■■■^ dy„, (3.15) 
h{y)---h„{y) 



s := ^deg/, -degp- (n + 1). (3.16) 

i=i 

Notice that r,hi,...,h„ e C[y]. 

Let [p] = [po ■ ■ ■ • : Pn] S LZj n V(Fi, . . . , F,;). We can define the residue of the form Q at 
point [p] as 

res(a[p]) :=res(c5,0,([^])), (3.17) 

where a; denotes form Q expressed in the chart (U,, 0,). 

The form Q is defined by homogeneous polynomials Fi, . . . , F^ and 

P(zo,Zi,...,z„):=Zo^^SPp(| I). (3.18) 

To underline the explicit dependence of O on F, and P we write symbolically Q = P/F. 
The following theorem is a special version of the global residue theorem. 

Theorem 3.7. Let V{F) := V(Fi, . . . , F„) be a finite set. Then for each polynomial P such that 

degP<f^degF,-(n + l), (3.19) 

! = 1 

we have 

X] res(P/F,[s]) =0. (3.20) 

WeV(F) 



21 



For the proof and the more detailed exposition we refer the reader to 
If G V(/) is an isolated but not simple point, then we cannot use formula {3.9} to 
calculate the residue of the form co at this point. In such a case we can apply a very nice 
method developed by Biernat in [5; ^] that reduces the calculation of multi-dimensional 
residue to a one dimensional case. We describe it shortly below. 
Let us consider the following analytic set 

A:={xeU\ f2{x) = ---=fn{x)=0}, (3.21) 

where LJ C C" is a neighbourhood of the origin. Set A is a sum of irreducible one di- 
mensional components A = AiLS ■ ■ ■ LI Am- Let t (Pj{t) G Aj, ^,(0) = 0, be an injective 
parametrisation of Ai. Then we define the following forms 



As it was shown in we have 



res(a;, 0) = ^ res(a;i, 0) . (3.23) 
1=1 



3.3 Proofs 



Let V G Cj-[x], where x = (xi, . . . , be a homogeneous polynomial of degree k > 2. We 
define n polynomials / = {fi, - ■ ■ ,fn) putting 

dV 

fi = ^- Xi, for / = 1, . . . ,n, (3.24) 

aXi 

i.e., f{x) = V'{x) — X. In our considerations the basic role plays the following n form 

CO = ^ , , , dxi A • ■ ■ A dxn, (3.25) 

/iW • ■■fn{x) 

where p is a polynomial. The polar loci of this form are points of V(/) := V(/i, . . .,f„). Let 
us notice the following facts. 

Proposition 3.1. Point G V(/) is a simple point and /'(O) = —E„. Thus we have 

res(a;,0) = (-l)XO). (3.26) 

Ifde V{f) and d ^ then 

1. point [d] G CP"^^ is a proper Darboux point ofV, i.e., [d] G D*{V), 

2. the Jacobi matrix f' (d) has eigenvalues Ai{d), A„^i{d), A„{d) =k — 2, 

3. ifdetfid) ^ 0, then 

- (t-2)A.(y-A„_.W 

4. points dj := dd G 'V(/), where e is a primitive {k — 2)-root of the unity, satisfy f [df) = 
f'{d),forj = 0,...,k-3. 
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An easy proof of the above proposition we left to the reader. 

Let us notice that fj^ = diV, so /, intersect at the infinity at a point ti 7^ iff V (ti) = 0. 
For such d, point [d] G CP"^^ is an improper Darboux point. Hence, in order to investigate 
the intersection of / at the infinity we extend the form oj into a rational form O on CP". 
We proceed exactly as in the previous subsection. Namely, let 

F,-(zo,zi,...,z„) ■.= zf^f'f, f^,...,^) =^-4-2^,, (3.28) 

" V^O Zq J dZi ^ 

be the homogenisation of for i = 1, . . . ,n, and 

P(zo,Zi,. . .,z„) := Zq^'^^p (—,...,— 

the homogenisation of p. The affine part of the projective algebraic set V(F) := V{F-[, . . . , f„) 
is homeomorphic to V(/). Notice that V(F) = D{V). Hence, the form Q = P/F on the 
chart {Uo,6o) is given by (|3.25|) . According to the previous subsection on the chart (Lfi,0i) 
it reads 

^ = , I , dyi A • ■ ■ A dyn, (3.29) 

where 

s := n(/:-l) -degp- (n + 1), (3.30) 

and polynomials hj are defined by l|2.58|) and (|2.59|) . and r{y) is given by (|3.14|l . 

In the proofs of Theorems I3.1[ 13.31 and 13.41 we use the same idea: we calculate residues 
of the form O = P/F with appropriately chosen P at all points of V(F) = T){y). 



3.3.1 Proof of Theorem [ 

At first we prove the following lemma. 

Lemma 3.1. If assumptions of Theorem 13.11 are fulfilled, then for a polynomial P G C[q], with 
degP = m < n{k — 2) — 1, satisfying 



P{x) = P(e'x) for all x G C" and i = 0,...,k- 3, (3.31) 
where e is a primitive (k — l)-root of unity, we have 

_m 

[<i]GV 



Proof. We take the form CI = P/F and we calculate residues of this form at all points of 
V(F) = S(y). In the affine part of CP", form Q is given by (13:25)) . 

Assumptions of Theorem 13.11 imply that V'{x) = only for x = 0. Thus, polynomials 
f\, ■ ■ ■ , fn do not intersect at infinity. So all points of V{F) lie in the affine part of CP". 

If d G V(/i, . . . ,/„), and d 7^ 0, then [d] is a proper Darboux point of V. By our 
assumption, [d] is a simple point of T>*{V). By Lemma l2!6l all points of the preimage 
7r^^([d]) are simple. This is why the inequality holds 

n-l 

det f'{d) = {k-2)Yl Md) ^ 0. 

Moreover, by Proposition 13.11 dj := dd G . . . ,/„), for = 0, — 3, where £ is 

a primitive {k — 2) -root of the unity, and f' {dj) = f\d) for = 0, . . . ,fc — 3, so all Jacobi 
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matrices /' {dj) have the same eigenvalues Ai (d), . . . , A„_i (d), A„ (d) = k — 2. Additionally, 
polynomial p satisfies restriction (|3.10|) of Theorem l3.6l and p{dj) = p{d) for j = 0, . . . ,k — 3. 

The above considerations show that all assumptions of Theorem l3.6l are satisfied. For a 
given d G V(/i, ...,/„), d 7^ 0, we have 

rres(a;,d,) = — 3^^^^^ (3.33) 

;r"o Ai(d)---A„_i(d)' ^ 

and for G V{fi, . . .,fn) we obtain 

res(a;,0) = (-l)"p(0), (3.34) 

as /'(O) = —En- Taking the sum of (|3.33|) over all Darboux points and adding (|3.34|) we 
obtain zero. □ 

The elements of matrix f {x) are polynomials of degree k — 2. Thus, we can take 

p{x) = {i:rf{x)-{k-l)y, with r G {0,...,n-l}, (3.35) 
and apply Lemma [3.11 For this choice of p{x) we have 

p{d) = Ti{A{d)y for d G (3.36) 

and 

p{0) = {-n-{k-2)y. (3.37) 

Thus, making the above choice for p{x) we obtain relations l|3.2|) in Theorem 13. II 

In order to prove that relations H3.2|l are valid we define polynomials Pi{x) which are 
coefficients of the characteristic polynomial of matrix f'{x), namely 

det(/(x) - AEn) = {-!)" f A" - pi{x)A"-' + ■■■ + {-l)y„{x)) . (3.38) 



Obviously degp, = i{k — 2), and Pi{£'x) = Piix) for all x G C", j = 0, . . . ,k — 3, and 
i = 0, . . . ,n. For d G T)*{V) we have Pi{d) = Ti(A(d)), where 

A{d) = (Ai(d),...,A„_i(d),fc-2), 

and 

pr{o) = i-iy 

Applying Lemma ISHI for p{x) = Pi{x) with < i < n — 1, we obtain the following relations 
Since 

T,(A(d)) = T,(A(d)) + (/c-2)T,_i(A(d)), z = 0,...,n-l, (3.40) 
where T_i(A((i)) = 0, we have the following system of linear equations 



T, + (fc-2)T,_i = To = {-ir-\ 



(3.41) 
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where r = 1, . . . ,n — 1, and 



Tr= L -^'W^- (3-42) 



Now using induction with respect to r and the relation 

n — i \ / n — i — 1\ / n — i — 1\ 
r-i + l) ~ \ r-i J ^ \r-i + l)' 

it is easy to show that 

T, = {-iy+"-^p^(^~^_~^^{k-iy, for r = 0....,n-l. (3.43) 
This shows that relations (|3.2p are satisfied and the proof is finished. 

3.3.2 Proof of Theorem [3j 

Our proof is based on the following lemma. 

Lemma 3.2. If assumptions of Theorem 13.31 are satisfied, then for a polynomial p G C[^], with 
degree deg p = m < n(k — 2) — 2 = n — 2, we have 

^ T ?Ar^)) = (3-44) 

Proof. We proceed in a similar way as in the proof of Theorem 13. II However, now polyno- 
mials f\,. . . ,fn can have intersections at the infinity, and because of this we have to apply 
Theorem 13.71 Thus, we have to calculate residues of the form Q = P/f at all points of 
V{¥) = T){V). Points of T){V) which lie in the affine part of CP" are given by V{f). Since, 
by assumption, k = 3, there is one to one correspondence between points of T)*{V) and 
V(tf) \ {0}. Thus for each [d] G Ti*[V), by assumption, [d] is a simple point of T){V), so by 
Lemma IZ6l 7r~^([tl]) is a simple point of Ti{V). As a result for each non-zero d G we 
have 

det/'(d) = Ai(d) ■ ■ ■ An-i{d) ^ 0, 

and 

Moreover, the local residue of the form a; at x = is 

res(a;,0) = (-1)"P(0). (3.46) 

Let us assume that [d] G V(F) lies in the hyperplane at the infinity. Without any loss 
of the generality we can assume that [d] = [0 : di : • ■ ■ : dn] with di = \. We know 
that [d] = [di \ ■ ■ ■ \ d„] is an improper Darboux point which is, by assumption, a simple 
point of Ti{V), and thus, by Proposition 12.41 [d] is a simple point of T>{V). It turns out 
that we can calculate the residue of the form Q using the local residue formula on the 
chart {Ui,9i), where form Q is given by (|3.29|) . Since coordinates of [d] on {Ui,9i) are 
a = (0,^2/ • • - fdn) = 9i{[d]), and s given by formula (|3.30|) is greater than zero, we have 

res(a;, a) = 0. 

The above shows that points of V{F) which lie at infinity do not enter in the total sum of 
the residues. □ 
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To prove the theorem we apply Lemma l3!2l taking p{x) = (Tr/'(x) — l)'^ and p{x) = 
Pr{x) with r = 0, . . . ,n — 2. Then repeating arguments used in the proof of Theorem 13.11 we 
show that relations i3.2} and (|3.3[) are satisfied. 



3.3.3 Proof of Theorem HJ 

The proof will follow directly from the following lemma. 

Lemma 3.3. If assumptions of Theorem 13.41 are satisfied, then for a polynomial p G C[q], with 
degree degp = m < (n — 1) (A: — 2) — 1, satisfying 

P{x) = P{e'x) for all x G C" and i = 0,...,k- 3, (3.47) 

where e is a primitive [k — 2) -root of unity, we have 

Proof. As in the proof of Lemma |3^ we have to show that residues of the form O = P/F at 
points [d] G V(f ) which lie on the hypersurface at the infinity, vanish. The difficulty that 
we meet here is that now these points, by Proposition I2.5[ are not simple and we cannot 
use the local residue formula. 

Let us assume that [d] G V(f ) n Hoo. Without any loss of generality we can assume 
that [d] = [0 : 1 : d2 : ■ ■ ■ : d„]. The coordinates of [d] on the chart {Ui,9i) are a = 
(0, d2,...,dn) = 9i{[d]). We have to calculate the residue of the form co given by (|3.29|l at 
point a. We know that [1 : ^2 : ■ ■ • : d„] is an improper Darboux point. By assumption it 
is weakly degenerated. To calculate res(a), a) we apply the Biernat formula, see the end of 
Section |3^ To this end we need to determine all the branches of the analytic set 

A:={yeUa \ h2{y) = ■■■= hn{y) =0}, (3.49) 

where Ua is a neighbourhood of a. Since a are the coordinates of an improper Darboux 
point that is minimally degenerated, by Proposition 12.51 we have 

rank ^J)^^' fa) = n - 1, 

d{yi,...,y„} 

and moreover 

a(/z2,...,//„) 

d{ij2,...,y„) 

Hence, by the implicit function theorem there is only one analytic branch passing by point 
a. We can always choose coordinates such that a = 0. Then, because 

(y) = yiVi'^' for 2 = 2,..., n, 

dyi 



and 



we can put 



dv 

— (0)=0, for i = 2,...,n, 
dyi 

t/l(^) = ^, and yi{t) = 0, for 2 < i < n, (3.50) 
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as the parametrisation of the unique branch of A. Then, the Biernat formula (|3.22|) , applied 
to the form co gives the differential form iv{t)dt, where 



w{t) 



r{y{t)y 



Calculations show that 
where 

As an effect we have 



h^{y{t))deih'{y{t))dt 
deih'{y{t)) = -(fc-2)rt(f). 



h{y{t)). 



(3.51) 
(3.52) 



h{t) = det 



(0) - t'^-%j 



2<i,j<n 



w{t) 



r{t)t 



s-(k-2) 



h{t) 



where r(f) = r(y(f)) is a polynomial and moreover 



h{0) = det 



dyjdyj 



(0) 



(3.53) 



(3.54) 



2<i,j<n 



by assumption that the considered point is weakly degenerated. Thus, ^(f) is regular at 
t = 0, ifi s > k — 2. Thus, ii m = deg p < {n — 1) {k — 2) — 1, then iu{t) is regular at f = 0. 
But this is the assumption of our lemma, so all the residues at points vanish at infinity. 
We calculate the residues at points in the affine part of V(f ) as in the proof of Lemma ISJl 

□ 

To prove the theorem we have only to check if we can apply Lemma 13.31 for p{x) = 
(Tr/(x) - {k- 2)y with r = 0,. . .,n - 2. But deg(Tr/'(A:) - {k - 2))' = r{k - 2), so 
deg p < (n - 1) (/: - 2) - 1 for r = 0, . . . , n - 2. 

3.3.4 Proof of Theorem [3l2l and Theorem [33] 

Let us recall that A„^j( denotes the set of all admissible (n — l)-tuples for a given k. We 
define the set of distinguished spectra J^j^, which appears in Theorem 13.21 and 13.51 as the 
set of admissible (Ai,. . .,An) G A^^. satisfying all relations (|3.2|) and (|3.3p. In particular, 
(Ai, . . . , An) G 3^]^, satisfies relation (|3.3|) with r = n — 2, i.e.. 



^ T„_2(A0 _ (fc-l)"-n(fc-2)-l 
If we denote (Ai, . . . , Ajv) = (Ai, . . . , Am), where M = {n-\)N, then (|3.55P reads 

We need the following technical lemma. 

Lemma 3.4. for k > 2, n > 2 and N > 1, the set of solutions of (|3.55|) in A^j^ is finite. 



(3.55) 



(3.56) 
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Proof. We rewrite relation (|3.56|) in the form 



M 

^X; = -c, where c> 0. (3.57) 

As it is easy to check, the admissible set C Q of values of X, for a given k > 2, has the 
following properties: X;c = U X^, where X^ is finite and its elements are negative; set 
X^ is infinite and its elements are positive. From (|3.57|) it follows that a certain number of 
Xj must be negative, and we have only a finite number of choices of these negative X, . For 
each of them we have to investigate equation 

P 

X,- = c', where c' > 0, (3.58) 

! = 1 

where X, G X,^, for z = 1, . . . , p < M. But the set X^ has an important property. Namely, 
its only accumulation point is zero, i.e., X^ \ X^ = {0}. Hence, applying Lemma B.l from 



|220 we obtain that the set of solutions (Xi, . . . , X„) G {Xt)P is finite. □ 



Now, we pass to the proofs of Theorems l3.2l and l3.5l If an integrable potential V satisfies 
assumptions of one of Theorems 13.11 1331 or 13.31 then 

L{V) := (A(d) I [d] G D*(y) ) G A^^^, where N = card'D*(y), (3.59) 

fulfils, among others, relation (|3.3|) with r = n — 2. By Lemma I3.4[ this relation is fulfilled 
only by a finite number of admissible elements. 

3.4 Discussion 

We know that, for a given k > 2, and n > 2, the set of distinguished spectra 3^^. is finite 
but we do not know how many elements it has. We show that this set is not empty. 

Proposition 3.2. For each k > 2, n > 2, and each N > 1, set 3^^. is not empty and it contains at 
least two elements. 

Proof. We will use the following convention. If a collection 3 G C^j^ contains ni copies of 
the element ci, ... , and ny copies of the element ey, then we write 

a = (ni X ei, . . . , nj x Cj) = {ni x e,- 1 z = 1, . . . ,;), 

instead of 

a = (ei,...,ei, ...,e/). 

ni times 

For given k and n we define 

tti := ( -1, . , -\ , k -2,..^,k- 2 ), for z = 0, . . . , n - 1, 

n — \ — i times ; times 

n—2—i times 

c, := ( -1, . .J , -\ , k -2,...,k- 2 , k + 1), for z = 0,...,n-2, 

times n—i—2 times 
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and 

3i:= (^{k-2y(^."^ xuil i = 0,...,n-iy (3.60) 
32 ■= [oLi X ai,^iXbi,jiX Ci\ / = 1, . . . , n - 2; f = 0, . . . , n - 2 ) , (3.61) 

where 

«, :=(*-2)"-'-^(;:i), ft 2)— '("-'), 

7, = (t-l)(t-2)"-'-^(";^). 

Now, it is only the matter of simple but lengthy calculations to check that 3i and 32 are 
elements of with N = D{n,k). □ 

One can wonder whether there exist potentials V with spectrum given by 3-i or 32- 
Below we give examples of potentials satisfying this condition. 
Let us take 

Vi = f^Viqi where . . . , y„) G C". (3.62) 

i=l 

For generic values of parameters this potential possesses D(n,k) Darboux points and 
one can check that (A(d) | [d] G ©(Vi) ) = Si- 
Potential 

m] /, _ -N 

^2=E2""' ■ )'^i4'^' + L<^rql where (cg c„) G C"-^ (3.63) 

(=0 V ^ / ;=3 

possesses also D{n,k) Darboux points (for generic values of parameters). Moreover now 
(A(d) I [d] G V{V2) ) = 32. 

It is obvious that one would like to have theorems more general than those proved in the 
previous subsection. Having relations of the form given in the theorems, we have to keep 
the assumption that all proper Darboux points are simple. Hence the only possibility is 
to weaken assumptions concerning improper Darboux points. However, due to theoretical 
difficulties, we were not able to find good generalisations of the theorems. For example, 
we tried to relax the second assumption of Theorem 13.41 and proceed in the same way as in 
its proof. However we meet a serious difficulty. We did not find a general enough method 
to distinguish and parametrise branches of the analytic set (|3.49|) . 

Examples suggest that it is possible to find the desired generalisation. 

Example 3.1 Potential 

V = cjiiaql + hq^q^,), with a ^ 0, and b ^ 0, (3.64) 

has five proper Darboux points [dj], i = 1,. . .5 and two improper Darboux points. All of them are 
simple, and 

, / 3a + b 3a-b\ , , / 2{3a-b)\ 
A(di) = [-^^-^ ) ^ Md2) = Aids) = ^-2, ^ j , 

AW = A(d5) = (-2,-^^^). 

Although V does not satisfy assumptions of Theorem 13.11 all relations (|3.2|) are fulfilled. Of 
course (|3.64|) satisfies assumptions of Theorem 13.31 
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Example 3.2 Potential 

V = qiiql + iqiqi + aql), 

has, for generic values of parameter a, four simple proper Darboux points and one improper Darboux 
point which has multiplicity 3. Although assumptions of Theorem \3.3\ are not satisfied, relations (|3.2|) 
and for r = 0,1 are fidfilled. But 

y (Ai(rf) +A2(rf))^ _ l + fl(28fl-9) 
l,]ehiv) Md)A2{d) - 2^^ • 

Example 3.3 Potential 

V = {q2- iqif{biql + b2q\ + bsql + fo4^?i'?2), 

has, for generic values of parameters, three simple proper Darboux points and a curve of improper 
Darboux points [1 : i : a], where a S C. None of our theorems can be applied for this potential. In 
fact, we have 

y 1 ^ 4bib2+2{bi + b2)b3-bl 

[d]hiv)Md)A2{d)- 2{b, + b2Y 

^ Ai(rf)+A2(d) ^ 2{bl + bl) + bl 
[i]k^(v) Md)A2{d) 2{b,+b2Y ' 

^ (Ai(rf)+A2(tf))^ ^ 4(bf + b^b2 + bl) + 2(bi + ^2)^3 + hi 
[,]h{v) Ai(d)A2(d) 2{b, + b2Y 

so none of the relations of the form l|3.2|) is satisfied. Nevertheless, we have one universal relation for 
the considered potential. Namely combining the above formulae we obtain the equality 

y {l+A^{d)+A2{d)f _^ 

[,]kiv) Md)A2{d) 



4 Three degrees of freedom and potential of third degree 

In this section we consider systems with three degrees of freedom with a homogeneous 
potential of the third degree. The aim is to show how to apply our theorems to distinguish 
integrable potentials. In fact, we consider here only generic potentials having a maximal 
number of simple Darboux points. Investigations of non-generic potentials will be pre- 
sented in a separate publication. 

The whole analysis can be divided roughly into two steps: 

1. determination of the set of distinguished spectra d^j^, 

2. reconstruction of the potential. Knowing a distinguished spectrum (Ai,. . .,A]v) G 
•^nk have to find all potentials V with N proper Darboux points [dj], such that 
A,' = A(di)forz = l,...,N. 
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4.1 Distinguished spectra 

For n = k = 3 we have N = D{n,k) = 7, so we look for seven admissible pairs Aj 
(A^'-*, Aj"*) G ^3,3, z = 1, . . . ,7, which satisfy the following three relations 



1 



7 

SaPa^"^' 

f AP+Af _ 



! = 1 



16. 



(4.1) 



Numbers Aj'^ + 1 belong to the set M3 which corresponds to items 1-6 in the table of the 
Morales-Ramis Theorem ll.2l i.e.. 



M3 



-p{p - 1) I peZ\\J 



+ 3p(p + l) 



p G Z 



u 



.1 + A(l + 5p)2 
24 50 ^ 



pez|u|-^ + A(2 + 5pf I pGzj. 



The first step in solving (|4.1|) is to find all solutions of the second relation which we write 
in the form 

14 

X]X/ = -4, (4.2) 

/=i 

where X; = 1 / A^ . Knowing M3 we easily determine the set X3 of admissible values of 
X/. It has, as it is easy to check, the following properties: 

1. X3 = X3 U and X3 is finite, 

2. if X G X3 , then x < -1, 

3. if X G X^, then x < 1. 

It follows that among X; satisfying (|4.2|) . there is a certain number of negative ones. More- 
over, at most nine of X/ are negative and we have a finite number of choices of them. 

Example 4.1 Assume that Xi,...,X7 are negative. Taking Xi,...,X7 G X3" we obtain 2426 
different choices satisfying Xi + ■ ■ ■ + X7 < —4. 



Thus we reduce the problem to finding solutions of 

c, where c > 0, p < 14, 



/-I 



(4.3) 



If (Xi,. . .,Xp) is a solution of the above equation, then we can assume that Xi < X2 < 
• • • < Xp. Hence Xp > c/p, and we have only a finite number of choices for Xp. For each 
of this choices we have to find the solution of equation (|4.3|) where we replace p by p — 1. 
Repeating this reasoning we end up with equation 



Xi + X2 = e. 



(4.4) 
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Table 4.1: The distinguished spectra for n = k = 3 

1. (3x (-l,-l),3x (-1,1),(1,1)), 

2. ((-1,-1), (-l'-0 ' (l'-0 '2 X (-1,4),2 X (1,4)), 

3. ((-1,-1), (-1'-^ ' (l'-0 '2 X (-1,14),2 X (1,14)), 

4. ((-1, -1), (-1,-0 , (l, -I) , (-1, , (l, I) , (-1, 14), (1, 14)) , 

^ ((4-D'(-^)'(-^)-(H-<-^0' 

7. (2x(-l,6),(-l,0,(-^.-|),(H.14),2x(14.39)), 
"!'"§) '^^ (-l'5)'3x (6,14)), 

10. ((-l,l),(-l,-^),2x(-|0,(l,H),2x(14,39)). 

SO, e > X2 > e/2, and we have only a finite number of choices for X2. This procedure 
proves that we have only a finite number of solutions of (|4.3|), however a serious practical 
problem appears because we do not know a priori the lower bound for Xi . We can order 
elements of into decreasing sequence {x,}. Then Xj decreases with i as i^^. For a small e 
interval [e/2, e) contains approximately [1/e] admissible elements X2. For each of them we 
have to check if Xi = e — X2 is admissible, i.e., if Xi G X^. Examples show that e can be of 
order 10^^^ and smaller. Very small values of e in (|4.4|) appear when the number of negative 
elements is smaller than seven. For such cases, even using a dedicated software rurming 
on the quickest accessible computers, we were unable to perform the computations up to 
the end in a reasonable period of the time (one month). This is a reason why we restrict 
our searches giving a priori lower bound for the smallest element Xi . 

A solution of (|4.2|) gives 14 numbers (Li, . . .,£14), where L, = 1/X,-. We divide them 

into all possible seven unordered pairs A, = (a['\ A2 ^) and check if these pairs satisfy the 
remaining two relations (|4.1|). 

Using the described algorithm we have found the distinguished spectra listed in Ta- 
ble O 



Remark 4.1 The described algorithm has an obvious generalisation for arbitrary k and n. 
4.2 Normalisation of the potential 

The C-linear space Cj; [q] of homogeneous polynomials of degree 3 in three variables has 
dimension 10. Starting from the general potential of degree three 

V =aiq\ + a2q\q2 + a^q\q^ + a^qiql + asql + aeqlq3 + ajql + agqiql + flg'?!'?! r^ 

(4.5) 

+ awqiq2q3, 
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we have to select its 'good' equivalent representative with the minimal number of free 
parameters. The group PO(3, C) is four dimensional, so a typical orbit of its action on 
C3 [q] has dimension 6. However, the problem is that there is no good parametrisation of 
these orbits. 

Here we consider only a generic case and this simplifies the analysis a lot. In our further 
considerations we will need the following fact. 

Proposition 4.1. Assume that a potential V G C3[^] has a proper non-isotropic Darboux point [d] 
such that matrix V"{d) is semi-simple. Then it is equivalent to 

V = aiq\ + a2q\q2 + fl3^?i'?3 + «4'?i'?2 + «5^?2 + (^e^h?, + ^^l- (4.6) 

Moreover, d = (0,0, 1) and matrix V"{d) has eigenvalues Ai = la^, A2 = 2a(, and A3 = 2. 

Proof. By Remark 12.11 we can assume that d = (0,0,1). This implies that in (|4.5|) we have 
ag = ag = and ay = 1/3. Vector d is an eigenvector of V"{d) with eigenvalue A3 = 
A: — 1 = 2. Hence V"{d) has the following form 

2fl3 «io 
«io 2fl6 
2 



V'{d) 



(4.7) 



and its remaining eigenvalues are Ai^2 = ^3 + ^6^ \/ ^10 + (^^3 ~ ^e)^- If 

flio 7^ ±i(fl3 - fle), (4.8) 

then we can make a rotation around d such that in the new frame flio = 0. Under this 
assumption the normalised potential reads 

V = aiq\ + a2qlq2 + a^qlq^ + a^qiql + asql + 06^72^3 + ^^3- (4-9) 

If condition (|4.8|) is not fulfilled, then V"{d) has a double eigenvalue Ai^2 = ^3 + ^e- In this 
case, because V"{d) is semi-simple, we have 

rank(y"(d) - (03 + ae)E3) < 2. (4.10) 

This condition gives 

(fl3+fl6-2)(fl3-fl6) =0. (4.11) 

If fl3 + flg — 2 = 0, then all eigenvalues of V"{d) are 2. But now, if 03 ^ 1, then rank(y"(d) — 
2£3) = 1. Hence, 03 = flg = 1, and so flio = 0. In the alternative subcase 03 = fle, so again 
flio = 0, and this finishes the proof. □ 

Notice that for the normalised potential (|4.6|) the proper Darboux point = [0 : : 1] 
lies in the line at infinity Hco. We called it the Darboux point at the infinity, although, as 
the following proposition shows, we can have more Darboux points in this line. 

Proposition 4.2. Potential (|4.6|) has at most 3 proper Darboux points in the line at infinity Hoa- if 
fl4 7^ 0, then it has exactly one proper Darboux point at infinity. 

Proof. If [q] is a proper Darboux point of V which lies in Hoo, then V\q) = q and q-[ = 0. 
For the potential (|4.6|) this gives the following system of equations 

a^q\ = Q, (?2(-l + 3fl5(?2 + 2^6(73) = 0, a(,q\- q^, + ql = 0. (4.12) 

It is easy to see that if a^ ^ 0, then the above system has only one non-zero solution 
corresponding to [rf] = [0 : : 1] . Thus, 04 = 0. Now, we have equations 

- 1 + 3fl5^?2 + 2fl6i73 = 0, a(,q\ -qo, + ql = Q, (4.13) 

that admit at most 2 solutions. □ 
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4.3 Reconstruction of potential 

At this point the problem is following. For each distinguished spectrum given in Table 14.11 
we have to find aU non-equivalent potentials having such a spectrum. This problem is dif- 
ficult. We have to determine the coefficients of the potential knowing only the eigenvalues 
of its Hessian matrix calculated at certain points with unknown coordinates. 

For each distinguished spectrum from Table 14.11 we can assume that the potential has 
the form (|4.6|) . So, we assign a pair (Ai,A2) from the spectrum to the Darboux point at 
infinity. By Proposition 14. 1[ we have Ai = la^ — 1 and A2 = 2fl6 — 1- In this way, we 
fix values of two unknown parameters fls and flg- Hence, the unknown parameters are 

We do not know the coordinates of the remaining Darboux points. So we have the 
following problem. Let us assume that we know eigenvalues Ai, A2 and A3 = /c — 2 of 
matrix f'{q) = V"{c]) — £3, at an unknown point q satisfying f{q) := V'{q) — q = 0. 
We have to derive restrictions on the coefficients of the potential. A direct approach to 
this problem is based on a simple idea. On the one hand, we know numerical values of 
symmetric polynomials s, = t, (Ai, A2, A3) for / = 1, 2, 3. On the other hand, we can express 
these values as coefficients pi{q) of the characteristic polynomial of matrix f'{q), namely 

det(/(^) - A£3) = -(A3 - \%{q) + Ap2(^) - p3{q)). 

As a result we obtain the following system of polynomial equations 

fi{q) = 0, pi{q) := Pi{q) - s,- = 0, / = 1,2,3, (4.14) 

where polynomials /, and p,- are considered as elements of the ring C [a] [q] . An elimination 
of variables q from the above system gives rise to a set of polynomial equations for the 
coefficients of potential V. We repeat the calculations for each Darboux point and then find 
solutions of all obtained equations. 

Remark 4.2 If the considered distinguished spectrum contains I < 7 different admissible pairs, 
then, because one of them is attached to the Darboux point at the infinity, we have at our disposal 
I — 1 pairs. Hence, ifl<7, we have to attach the same pair to two or more Darboux points. As an 
effect, we have to apply the algorithm only for I — 1 Darboux points with different attached pairs. 

The elimination step in the described procedure consists of determination of the so 
called elimination ideal of a given set of polynomials, see It can be performed only 

with the help of a computer algebra system. However, in practice, the known algorithms 
work effectively only in cases with a relatively small number of variables. Hence, in order to 
have a chance to obtain the desired result, we have to minimalise the number of unknowns. 
In the considered problem we can do this introducing the affine coordinates xi and X2 for 
the unknown Darboux point. They satisfy the following polynomial equations 

§■1(^1,^2) =0 and ^2(^1,^2) =0, (4.15) 

where 

Si-=^.-^iSQ' for 2 = 1' 2, 

and 

dv , , ^^.^ , 

gQ=]w- Xi— - 3^2^^, V[Xx,X2) = V{l,Xi,X2), 

see Lemma 12.11 and its proof. Moreover, we take values of two symmetric polynomials 
Si = Ti(Ai, A2) = Ai + A2 and S2 = T2(Ai, A2) = A1A2. Then, we have to express Si and 
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§2 as functions of xi and X2. To this end we recall that the eigenvalues of f'{(]) coincide 
with the eigenvalues of Jacobi matrix f{x) defined by (|2.39|) . Performing simple calculations 
like in the proof of Lemma [2.51 we find that 

^oS.^('S^'S) and ,§S.= f|l|§-|!l|!lV (4.16) 



dxi dx2 J \dx1dx2 dx2dxi 

Now, eliminating Xi and X2 from equations (|4.15|) and (|4.16|) we obtain a set of equations 
for unknown coefficients of the potential. 

Remark 4.3 The drawback of using the affine coordinates is following. We do not know in advance 
how many points are located at the infinity and what are their admissible spectra. Thus we have 
to check all the possibilities. One can perform another normalisation of the potential with the as- 
sumption that all Darboux points are located in the affine part od CF-^. However this causes other 
computational problems. 

Using the described algorithm it is always possible to derive a set of polynomial equa- 
tions which must be satisfied by the coefficients of the potential we are looking for. How- 
ever, generally the equations are very complicated, and it is impossible to find their explicit 
solutions even with the help of the strongest tools of the computer algebra. In order to 
overcome this problem we are forced to apply a different strategy. 

In fact we do not want to know the explicit form of all potentials with a given distin- 
guished spectrum, only those which are integrable have to be selected. Our potential satis- 
fies a priori all the necessary integrability conditions of Theorem ll.2[ However, the stronger 
integrability conditions were formulated by Morales, Ramis, and Simo, see [23; 24; 28]. 

Theorem 4.1 (Morales, Ramis, Simo, 2006). If a complex Hamiltonian system is integrable in 
the Liouville sense in a neighbourhood of a particular non-equilibrium solution, then the identity 
component of the differential Galois group Sm of m-th order variational equations along this 
solution is Abelian, for an arbitrary m G N. 

Generally it is difficult to apply the above theorem, because higher order variational 
equations have a big dimension. However, if the first variational equations are a direct 
product of Lame equations of the form 

d^Y 

= [n{n + l)p{t) + B]x, neZ, B G C, 

then we can use the following local criterion. If Si is not a finite group and a local solution 
around i = of m-th order variational equations contains a logarithmic term, then S^ is 
not Abelian, see Appendix A in the recent survey of Morales and Ramis [29]. 

For details how to apply this criterion for homogeneous potentials with two degrees of 
freedom see [2ll:l22l ]. In the case of homogeneous Hamiltonian systems with three degrees 



of freedom the strategy is completely similar to this applied in 11211: 12211 . We remark only 
that the criterion is applicable for a solution generated by such a Darboux point [d] that its 
eigenvalues A, of V"{d) are admissible and all are non-negative integers. In all cases, when 
we apply the higher order variational equation, we take a particular solution generated by 
the Darboux point at the infinity. 

4.3.1 Families 1-4 

Let us take as the distinguished spectrum the first line in Table 14.11 From Lemma 12.51 it 
follows that a potential with such a spectrum admits a proper and non-isotropic Darboux 
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point. We can assume that this point is [<i] = [0 : : 1] and assign to it the admissible 
pair (Ai, A2) = (—1, — !)• Moreover, by Remark ll.3[ we can assume that matrix V"{d) is 
semi-simple. As an effect, we consider potential of the form (|4.6|) . By Proposition 14. 1[ we 
have Ai = la^ — 1 and A2 = laj, — 1, so = Uf, = 0. But then in the potential (|4.6|) variable 
q3 is separated from qi and q2- This fact, and the assumption that the potential has seven 
Darboux points implies that one can make a rotation in the plane {^1,^2) and achieve that 
fl4 = 0. Then the potential simplifies to the form 

V = ami + mlqi + «5^?2 + 3^3' (4-17) 
and in the affine coordinates we have 

gl = (3fl5 - 2fl2)x? - 3fliX2 + fl2/ g2 = ^2- 3^2(301 + lUiXi) , 



3fli + 2fl2^i- 



(4.18) 



By the Bezout theorem, the number of solutions of =§^2 = is not bigger than four, so 
there is at most four Darboux points in the affine part of CP^. Since by Proposition I4.2[ at 
most three Darboux points can lie at infinity, and we assumed that the potential has seven 
Darboux points, four Darboux points must lie in the affine part of CP^ and three at infinity. 
This fact implies that fls 7^ and 3fl5 — 2fl2 7^ 0. For each Darboux point in the affine part of 
CP2 we attach the admissible pair (Ai, A2) G {(-1, -1), (-1, 1), (1, 1)}. Next we calculate 
symmetric functions Si = Si(Ai, A2), S2 = S2(Ai, A2), and then, using a computer algebra 
system, perform an elimination of Xi and X2 from the corresponding equations (|4.15|) and 
(|4.16|) . This elimination gives only one restriction on parameters fl„ namely fl2 = 0. But 
then, from (|4.18|l we see that = fli . By assumption all Darboux points are proper, and 
we have necessarily ai ^ 0. The affine coordinates of four Darboux points are following 

(^1,3:2) = (0,0), {xi,X2) = (0,3fli), {xi,X2) = ' (^l'^2) = ^^,3fli^ . 

We can check directly that the non-trivial eigenvalues (Ai, A2) are as we required ( — 1,-1), 
(—1,1), ( — 1,1) and (1,1), respectively. 

Assumption As 7^ guarantees that three Darboux points are localised at infinity, and 
their coordinates are following 

{(]\, (]2, qs) = (0, 0, 1), {qi, q2, qs) = (^0, 0^ , ((?i, (?2, ^3) = (^0, 1^ . 

The non-trivial eigenvalues (Ai,A2) at these points are (-1,-1), (-1,-1), and (—1,1), 
respectively. 

In this way we showed that if a potential admits seven proper Darboux points with 
a distinguished spectrum given by the first line in Table I4.1[ then it is equivalent to the 
separable potential 

Vi = aiql + asql + -ql, aiUs + 0, (4.19) 
with obvious additional first integrals 

11 1 

^1 = + 3^3' h = + fls'?!- 

An analysis for the next three distinguished spectra is completely analogous to that 
described above. For the spectrum given by the second line in Table |4T] we obtained one 
integrable potential 

^2 = Y^lq2 + asql + pi fls 7^ 0, (4.20) 
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with first integrals 

1 1 

The third spectrum gives rise to another integrable potential 

^3 = ^(ihl + «5'?2 + «5 ^ 0, (4.21) 

with additional first integrals 
1 1 

^1 = 2^3 + -^4' h = 128p\ - a5ql{32pip2qi - 96plq2 + asqliql + 6ql)). 
The fourth spectrum generates also an integrable potential 

^4 = asiiVSql + 9(??(?2 + I8ql) + fls 0, (4.22) 

which possesses two additional first integrals 

1 2 1 3 

h = 3p\ + 2i\/3p\p2 - 54,a^pip2q\{qi + i\/3(^2) + \8a^p\qi{i\/?>q\ + 6^?l^/2 + 6i\/3q\) 
+ 9ia5ql{2V3pl + 3fl5(5i(?f + 6v/3^/?(?2 + 9i(?i(?^ + uVSqD). 

Let us note that potentials Vi,V2,V2, and V4 have the form Vi{qi,q2,q3) = Vi{qi,q2) + 
ql/3, where Vi{qi,q2) are integrable two dimensional homogeneous potentials. They have 
appeared in the famous Hietarinta table in [15] and later were re-obtained in the differential 
Galois framework in [21] where it was shown that they are the only integrable potentials 
with the maximal number of Darboux points for k = 3. 

The reconstruction of the potentials corresponding to the next families is more compli- 
cated by reason of the absence of the admissible pair ( — 1,-1). Still one can easily decode 
from the pair of As at the infinity two coefficients of the potential but, separation of the part 
of V depending on {qi,q2) and qj no more appears, thus we cannot put ^4 = 0. Thus all 
six Darboux points rest at the affine part of CP and the level of computational difficulties 
increases very quickly. 

4.3.2 Families 5-7 

For the fifth distinguished spectrum in Table |4T] we assign the pair (Ai, A2) = (—1,4) to 
the Darboux point at the infinity. Then we obtain immediately two coefficients fla = and 
fl6 = 5/2 of the potential l|4.6p. Thus, we will analyse the higher order variational equations 
for the potential 

5 1 

V = aiql + a2qlq2 + fl4(?i(?2 + ^sql + 2'?2'?3 + ^ql- (4-23) 

A local solution around i = of the second order variational equations does not have log- 
arithmic terms. Such terms appear in the solutions of the third order variational equations. 
Comparing them to zero we obtain the following equations 

a4(35fl2 - ISfls) = 0, -147 + 28fl| - 27fl§ = 0, fl2(9fli - ^4) = 0, 
84fl| + 168flifl4 - 56fl| - 27fl2fl5 = 0, 
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with solutions 



fl2 = 0, fl4 = 0, «5 = ±|- (4.24) 
fli=0, fl2 = ± J, fl5 = ±y, fl4 = 0, (4.25) 



i / 21 , / 14 9i / 21 , 35i / 14 



^^^ = -^269' ^^^ = ^^^269' '^ = -2V 269' '^ = ^XV 269' ^^"^^^ 



21 , / 14 9i / 21 , 35i / 14 



" " °' ''^ " «5 = 0. (4.28) 

The first solution (|4.24|) gives the integrable potential 

7i 5 1 

which is equivalent to potential V2 defined by (|4.20|) . 
Solution (|4.25|) generates the integrable potential 

V5 = jqlqi + jql + ^qhs + ^ql 

with first integrals 

h = 48pi(ip2 + P3)qi - (3^?? + 2{q2 - 2iq3f){3ql + 2q2{5q2 - Aiqs)) - 48p?(i^?2 + qs) 
+ 32{p2 - 2ip3){p2q3 - psqi), 

h = 24pf - 8p?(12(p2 - ip3)(P2 - 2ip3) - 9iqlq2 + 2{2q2 - iq3){q2 - 2i^j3) (5i(?2 + 4^?3)) 

+ 3ql{3q\ + 16(p2 + ips) (P3'?2 - P2'?3) + 6ql{3ql - 4i(?2(?3 + ql) 

+ 4.q2{q2 - 2iq3){7ql - 4i(?2(?3 + 2ql)) + 2Ap^qi{p3{-3ql - 22ql + 20i(j2'?3) 

- ip2{3ql + 6^/2 + 4i(?2'?3 + 8ql)). 

Solutions (|4.26|) - (|4.28|l give potentials 



- / 21 3 / 14 2 9i / 21 2 35i / 14 3 5 2 1 3 
^ = "2 V 269^1 ^ '^V 269^1^^ " 2 V 269^^^?^ ^ 269^^ + 2^?^^^ + 3^^'' 



i / 21 3 / 14 2 9i / 21 2 35i / 14 3 62 1 3 
^ = 2 V 269 ^1 ^ V 269 ^^^^ + 2 V 269^^^^ ^ ^ V 269 + 2^^^^ + 3^^, 



, 1 /7 3 , V21 2 5 2 1 3 
2 V 3'^^ ^'?i'?2 + 2'?2'?3 + g'Js- 

All of them are not integrable because the respective solutions of the fifth order variational 
equations have logarithmic terms. 

For the sixth distinguished spectrum we assign pair (Ai, A2) = (4,21) to the Darboux 
point at the infinity and we apply the same procedure as in the previous case. Solving 
conditions obtained from an analysis of solutions of the third order variational equations, 
we obtained only one integrable potential 

Ve =364Vl7(j? + 2S35W\7qlq2 + 1560v/l7^?irj^ + 6552i\^ql + 4335rjf (?3 

+ l9Q7iqlq3 + 578^?^, ^^'^^^ 
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with one first integral of degree two in the momenta 

h = MVrJflqi + 8p2P3(18i(?i + qi) - 2pl{7VV7qi + 4r?3) - p\{54i^q2 + 2(>q^) 
+ 2pi{p2{27iVl7qi + 7VUq2 + 9iq3) + P3(13(?i - 81i(?2 - IvV^Vqs)) 

- 51(-2568^j| + 36i^j^(695i/i + b2\fY7q3) - 36iqiq2{151ql + lAlVlVq^qs + 221(?§) 

+ 2^?i(31351rjf + 5386\/T7(ji(?3 + 2924(?§) + qi{12205ql + 166S\fV7qlq3 + 3978qiql + 1156^17^?^)), 

and the second more complicated of degree four in the momenta is given in Appendix. 

For the seventh distinguished spectrum in Table |4T] we assign pair (Ai, A2) = (—1,6) 
to the Darboux point at the infinity. The conditions that there are no logarithmic terms in 
solutions of variational equations of orders 3, 4 and 7 yield one integrable potential 

V7 = U\/7ql + 2mVUqlq2 + 330V7qiql + 9351^14(^2 + 5087qlq3 + 294ql 
One its first integral is of degree four with respect the momenta 
h = I9p\ + 47p| - UiVUplp3 + 16p?(8i\/2p2 + 7V7p3) + 294p2P3(- 4361^2^^3 
+ 6qiq2{-l95iV2q2 + 56v^<?3) + 6(??(17'?2 - 77i^q3) + ql{827q2 - 399iv/l4r?3)) 
+ 252pl{l2\/7q\ + q\{145W\Aq2 - 9Sq3) + 2V7q^{185ql - 98ql) + q2{6mVTiql 
+ 2303(?2'?3 - ^9W\Aql)) + \8pl{-llpl + UiVUp2P3 + 7(7p| - 136V7ql + A80iVUqlq2 

- 2756V7qiql - 27iVUql + 62721 v/2(?i(?2^?3 - 9898(/i(?3 + 784V7qiql + 34301 v/l4(?2(?3 

+ 4116^?!)) + pi{-65iV2pl - U7V7plp3 - 147^3(1216^?? + 24(??(-213iv/2<?2 + 28V7q3) 
+ 96qiq2{172q2 - mVTiqs) + q2{^l89iV2ql + 596^^/7q2q3 - 617AiV2ql)) + 21^2(141 v/2p§ 

- 16641 Vl4(jf + 6607 V7ql - 28665i\/2qlq3 - 3822\^7q2q3 - 617^iV2ql 

+ 384:ql{AV7q2 - 491^2(^3) + 6qi{-1555iVUql + 392(?2(?3 - 1666iyT4(?§))) + 294(27v/7p|(8(j? 
+ 22iV2qlq2 - 24q^ql + 9iV2ql) + 112{434qf + 96ql{-7iV2q2 + 12V7q3) + 6q\{ll9ql 

- 32^i\/Uq2q3 + 924(?§) + 12qiql{20UiV2ql + 1782V7qlq3 - 26251 v/2(?2(?3 - 50AV7ql) 

+ ql{2233i^ql + 5931v/7(?^i?3 - 143641 v/2(?2^?3 + 1134v/7rj|) - 8ql{6986ql - 3492i^qlq3 

- 77^9q2q3 + ^78Wliql) - 3q\q2{3626ql + A8iVUqlq3 - 3465(?2(?3 + 1386iVUq3))) , 
and the second of degree six with respect to the momenta is given in Appendix. 

4.3.3 Families 8-9 

For the eighth distinguished spectrum from Table WTll we assign pair (Ai, A2) = (6, 14) to 
the Darboux point at the infinity. The conditions that there are no logarithmic terms in 
solutions of variational equations of order four are following 

fl4 (-7291848638025 - 195171588288fli - 7947435510flifl4 + 20220446680fl| 
+ 315697519440fl2fl5 - 183025477350^^) = 0, 

fl4 (411968789569575 + 2262392768724fl^ + 1291180817430flifl4 + 339419628940fl| 
+ 1213154206380fl2«5 +4209585979050fli) = 0, 
2013340175fli + 114593850flf + 9791320flifli - 6923862fl|fl4 = 0, 

729732302775fli - 633757824fli + 43836570855^4 + 8866402020fl?fl4 + 4915882928fl^fl4 
+ 1092887040flifl| + 217032200fl| + 16783940880flifl2«5 - 5584905348fl2fl4fl5 = 0. 

(4.30) 
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The above system has very complicated solutions. In order to discard those which do not 
satisfy our assumptions about the number of Darboux points and their spectra, we decided 
to proceed as follows. At first we assume that there are six Darboux points in the affine 
part of CP'^, i.e., by Proposition I4.2[ we assume that 04 ^ 0. For three of them we assign 
(— |,— I), (—1,0 and (6,14), as the respective admissible pairs (Ai,A2). For each of 
these points we considered system of equations (|4.15|) and (|4.16|) . Then we calculated the 
elimination ideal with respect to variables [xXfX'i) and its Groebner basis. Next we join all 
these three bases and polynomials given by the left hand sides of equations (|4.30|) . The 
Groebner basis of all these polynomials is very simple. Equating to zero all polynomials 
from this basis we obtained equations 

al + 28431/320fl2 + 729/64fl5 =0, a\- Vlbl\V7a\ - 600/13 = 0, 

, (4.31) 
fl2fl5 - 145/117fl| - 540/13 = 0, = fl4 = 0. 

So we obtained a contradiction with our assumption that 04 = 0. Hence, only four points 
are in the affine part of CP^. Assuming that to three of them the attached admissible pairs 
are as above, we arrive to the same equations (|4.31|) . Their solutions give rise to only one 
integrable potential 

7 7 5i-\/3 7 9i-\/3 15 , 1 ^ 

^8 = 2^l'?3 ^^l'?2 ^^?2 + Y^2^3 + 3^?3' 

that admits the following first integrals 

Ix = l<oi\f?>-p\ + 3pi(6iA/3p2 - 48p2P3 - 32i\/3p3 + T7c^^ + 69iV3qlq3 + 15(?2(13i?f - 8ql) 
+ iV3q3{91ql + 16ql)) + 3qi{p2{-91ql + 45ql + 42iV3q2q3 + 71ql) - ip3{65\^q\ 
- 99\/?>ql - 480i(?2(?3 + lll\/?>ql)), 

h = 8619pt + 459pt + 90001 V3p^p3 + 21519p^ + mmq\ + 78p?(141p^ + 220iv/3p2P3 + 41p3 

- 369iv/3(?^ + 15A7qlq3 + 1595qlq3 + ll^l + 85iv/3^j2(-13(?f + ^l)) - 18pi(2649p§ 

- 791i\/?>ql - 6357qlq3 - 765qlq3 + 86ql - 5iV3q2{299ql + WOql)) + 20p2P3(- 15681 VSp^ 
+ 6075ql + 9675iV3qlq3 - 189(?2(65i?f + lAql) + 21iv/3(?3(65(?f + lUql)) + 6p| (215731 V3(?| 
+ lWllqlq3 - 104085qlq3 + 4782ql - 5i v/3(?2 (3939^?? + 7840ql)) + 780pi(?i (p3(13(?5 + 693(?| 
+ 262iA/3(?2'?3 - 56ql) + 3ip2(13v/3rjf - 35^/3(^1 + 32i(?2'j3 - ^V^ql)) + 12(28674^jf 

+ 251101 v/3^?^(j3 + 30758qfql + 79\7q{ql + 797ql + ^05q\{-l3q\ + 3Qql) - \5qlql{-650{)q\ 

+ 381ql) + ^iV3qlq3{-llS83q\ + 5959(?§) + m^q2q3{-ll97q\ + \l\8q\ql + 392^?^)). 

It does not end the analysis, because it is possible that to four Darboux points in the affine 
part of CP^ only two different admissible pairs are attached. To analyse all these cases we 
proceed as follows. We know that two additional Darboux points are in line at the infinity. 

Their coordinates and spectra (A^'^,A2'') can be found explicitly. Since now = 0, they 
can depend only on three parameters fli, a2, and a^. In fact they depend only on a2, and a^. 
So (Aj^'', Aj^^) and {A^^\ A^^'^) are not independent, and the following relations hold 

(8 + 15AW)Af) = 14(Af)-A«), 

(8 + 15Af = -14(a5'^ - A^')), (4.32) 
(14 + 15Aj^Vr^ = -14Af). 
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Now, it is easy to check that from the considered distinguished spectrum there is only one 
choice of two pairs satisfying the above relations, namely 



and this gives 



fl2 = ±5i— - and fls = ±9i— -. 

2 2 

For both sign choices equations (|4.30|l imply that fli = 0. In this way, we recover a potential 
equivalent to Vq. 

For the ninth distinguished spectrum from Table 14.11 we assign the pair (Ai,A2) = 
(14, 39) to the Darboux point at the infinity, and then we proceed in the same way as in the 
previous case. Logarithmic terms appear only in solutions of the seventh order variational 
equations. The coefficients of these terms are four very complicated polynomials depend- 
ing on a-[,a2,a4,a5. We join these polynomials and all eliminating ideals obtained for five 
Darboux points in the affine part of CP^. The Groebner basis of this ideal is very simple 
and it gives the following equations 

2 4609248 ^ 2 1395372 „ 9 621 ^ 

This system yields the following integrable potential 

V9 = 27W3990ql + 3726VT5qlq2 - 456i V3990(?i(?^ - miVlBql - l\15q\q^ 
- 3mQqlq2, - SQql, 

with very complicated first integrals of degree 4 and 6 in the momenta which are given in 
Appendix. 

Potential V9 was obtained under assumption that six Darboux points are localised in 
the affine part of CP^. The cases when a smaller number of points is in the affine part of 
CP^ we analyse in a similar way as for the eighth distinguished spectrum. It appears that 
such cases for an integrable potential are impossible. 



4.3.4 Family 10 

For the distinguished spectrum from Table |4T] we assign the pair (Ai,A2) = (—1,1) to 
the Darboux point at the infinity. Doing this we additionally assume that the admissible 
pair ( — 1,1) is attached to a non-isotropic Darboux point. The case corresponding to the 
isotropic Darboux point will be considered separately. The potential has the form 

V = aiql + a2qlqi + a^qiql + asql + qlq^ + -q^. 

The solutions of the second and fourth order variational equations do not have logarithmic 
terms iff fl2 = and 

fl4(-20 - 144flifl4 + 160fl| - 45fl§) = 0, fl4(-100 - 90flifl4 + 128al - 225fl§) = 0. 
The above equations have the following solutions 

{fl4 = 0}, jfli = ^fl4,«5 = ±^\/-25 + 8fl|l . 
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The first solution yields the potential 

V = aiql + asql + qlq3 + -q^, 

that is equivalent to (|4.19|) . 

Let us take the second solution with the sign +. Then fifth order variational equations 
yield 

5 

fl4 = or fl4 = ± — —. 

If fl4 = 0, then the corresponding potential has only one Darboux point, so we discard this 
case. The second solution yields the following integrable potential 

_ 4.V2ql bqiql , 1 3 
Vio - + ^ + ^2^3 + 3^3. 

with first integrals of degree 4 and 6 in the momenta 

h = llpi - 27(?2 - IS^K'?? - 4v/2(?i<?3 + Iql) + 4(6pf - 3pi + IS^Jlql - Iql) (3pi + Iql) 
+ llql{?,pl{Vlqi - 4(?3) + llvm{qi + Vlq^) - 2ql{12ql + Vlqiqs + 2ql)) 

- 12p2<?2(2p3(16^?? + 3ql + 8V2qiq3 - Aql) + sVlpiiqj + 4rj§)) 

- 12pl{2p3{2V2pi + P3) - 4((J2 - q3)q3{q2 + qs) - V2q^{5ql + 8ql)), 

h = 81t?|(2V2^?i + qs) + 216p2P3'?2(v^'?i + ^qs) + SAq^iipl - ^pl + 4^2^?? - 24.q\qo, 

- 6V2qiql) + 384p2P3^??^?2(3pi + 8V2ql + 8qlq3 - 2V2qiql) - 72p\{3pl + 2ql) 

+ 144p2P3^?2(P2 + 8<??(2v^^?i + 3q3)) + lA^p\{V2plp3 + 3V2p2q2ql - ^p^qliqi + v^^?3)) 

- 32(p^ + 12p4^?^?3 + I2plql{^pl + ^q^ql) + 32q\{3pl + 2ql)) - l2p\{Ap\ 

- (^pmqiildql + 9ql + 8V2qiq3 - Aql) + 9q\{2q\ + AV2qiq3 + ql) + 32V2ql{3pl + 2ql) 
+ Uplipl + Aqlq3 - V2q,{ql - 2ql)) + 6ql{9V2plqi + 2^?i(-6(j? + 2V2q,q3 + ql))) 

- lUq^{pl{7ql + 5V2qiq3 + 2ql) + 3(??(3p§ - 2q3{-2ql + 2V2q^q3 + ql))) 
-^8ql{pi{5V2q^+Aq3)+Aplql{8ql + 2V2q^q3+3ql) + 8ql{9plqi + ql{-6V2ql 

+ 4(?i(?3 + V2ql))) + 6pi{16V2pjp3 + 16plp3{8ql - bqiql + 3V2qlq3) + ^plq2{-l6V2ql 

+ 32(^i(?3 + v/2(3(?^ + ^D) + 3p3ql{9V2qj - 64:ql{V2qi + 2(j3) - 12ql{2V2q\ + 8(?i(?3 

+ V2^i)) + 12p2qi{-3\/2plql + 9tjii/| + 32(?f(j§ + ^l^l + 6i/i(?^ + V2q3))). 

Let us consider the case when pair (Ai,A2) = (—1,1) is attached to an isotropic Dar- 
boux point. In this case we assign pair (Ai,A2) = (14,39) to the Darboux point at the 
infinity and this gives ^3 = 15/2 and flg = 20. Let [d] be the isotropic Darboux point 
corresponding to admissible pair ( — 1, 1). We have two possibilities: 

1. either ^3 7^ 0, and then using a rotation around the third axis one can armihilate one 
component, e.g. rfi = and d = (0, s,is), where s 7^ 0, or 

2. £^3 = and d = (s,is,0), where s 7^ 0. 

In the first case the condition that d = (0,s,is) is a Darboux point yields 

59. i 

fl4 = 0, fl5 = -yl, ^^19- 
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Then the eigenvalues of V"{d) are ((2ifl2 — 15)/19,2,2) and V"{d) is not semi-simple. Thus 
by Remark [L3] this potential is not integrable. 

In the second case condition that d = (s,is, 0) satisfies V'{d) = d leads to s = but it 
was assumed that s 0. 

5 Final remarks 

Presented calculations show that, formulated at the beginning, classification programme 
for homogeneous potentials is very ambitious and meets a lot of non-trivial theoretical, as 
well as computational difficulties. Natural question arises about possible simplifications 
and improvements allowing to dream about effective calculations for n = 3 and > 3, as 
well as for n > 3. One possibility is to use a full power of the higher order variational equa- 
tions techniques. Here we used only the local criterion of the existence of logarithmic terms 
in solutions of higher order variational equations, which implies that their differential Ga- 
lois groups are not Abelian. In fact it is the only known way to use the higher variational 
equations. However, it seems that, for systems with homogeneous potentials, stronger ob- 
stacles for the integrability can be deduced from an analysis of global differential Galois 
groups of these equations. 

Another question is if we have already detected all possible relations between the eigen- 
values of Hessian V"{d) at different Darboux points. In fvl] Guillot mentioned that other 
relations exist, but their form is completely not obvious. If such additional relations exist, 
then they give new obstacles for the integrability and they limit the number of elements of 
the distinguished spectra. Thus the answer to this question is of great practical ii np ortance. 

Let us mention also that relations (|3.2|) and (|3.3|) were already presented in |33[|. How- 
ever in [33] only the cases with maximal number of proper Darboux points were analysed, 
corresponding here to Theorem 13.11 Furthermore, the proof of this result was based on 
results of Guillot [12] obtained with the help of highly advanced techniques [4] for which 
there are no obvious ways for generalisations. In present paper Theorem 13.11 was obtained 
ab initio and its proof is completely elementary. Additionally, using the introduced here 
approach of multidimensional residue calculus, we were able to analyse also some non- 
generic cases. 

The first integrals given for the presented integrable potentials were found by means 
of the direct method [16]. We have postulated the form of a first integral as a weight- 
homogeneous polynomial function with respect to the coordinates and the momenta. We 
have started from the lowest possible weight degree equal to 2, and increased this degree 
up to finding the first integral. Thus the presented first integrals are of the lowest possible 
degrees. All of them are irreducible polynomials. This follows from the results contained 
in [20]. First integrals for some potentials are very complicated functions. We only mention 
that theses expressions have already been simplified. 

Another problem is that some integrable potentials possess imaginary coefficients. In 
all cases it is possible to convert them into real potentials, but together with a simultaneous 
change of the signature of the kinetic energy form. There is an open problem how to find 
their equivalent real form, or how to prove that such a form does not exist. 

Also the problems how to improve the calculations of distinguished spectra that are 
the starting point to look for new integrable potentials and their reconstruction are worth 
considering. Our calculations described in the last section show that the one distinguished 
spectrum gives rise to a finite number of potentials. However, in some cases (the first 
four distinguished spectra) we obtained families of integrable potentials depending on 
parameters. Let us notice, however, that this occurs only in cases when potentials separate. 
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i.e., if they are sums of two integrable potentials depending on different coordinates. 
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Appendix 

For potential 

Ve =364:\/V7ql + l?,3bWV7q\q2 + 1560\/V7qiql + 6552iVV7ql + 4335qlq3 + l9Q74.qlq^ + 57Sql, 
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the second first integral is following 

h = 20ipt(2^(?2 + 5i<?3) + 120p^<?3 + 24iplp3{29qi + 5i<?2 - uVUq^) + Apl{-2ip2{5V^qi 

- T\\fV7q2 - ?,\c\^) + pz(2bqx + 247i(j2 + Vd'i\fV7qz)) + pi(682193<?^ - 572220(?| 

+ 8p§(134v/T7rji + 42i\/l7(?2 - 713(?3) + 537856 v/T7(?f(?3 - 5080620(??(?§ + 1895840 v/T7(?i(j^ 

- 1228828(?| + 432i(?^(8279(?i + 2566\/T7(?3) + 180<?i (24259^?? + 1136\/T7<?i(?3 + 50286(?i) 
+ 72iq2{17153ql - 9978VV7qlq3 + 39882^1 ^§ - 8092VV7ql)) + p?(590325(jt - 11864844^| 
+ 215380V 17qlq3 - 1647300qlql + 714408 Vl7(ji«?| + 1268132^^ + pl{56VT7qi + 96i-s/V7q2 
+ 76(^3) - 4pl{lQ3VV7qx + {2m)VV7q2 + 1600(j3) + 72i<7| (60775(ji + 28124^/17^3) 

- 12ql{129965ql - 18626VV7qiq3 - 22542ql) + 36i(j2 (22525^? - 22680Vl7^?ij3 + 310386(?i<7^ 

- 50864Vl7(j^) + 4p2P3(439i?2 - 3i(103^i + 35^(j3))) + 4p2P3 {2709i VV7qt + 202986i\/l7^| 
+ 10(?f (20909^17(^2 - 158661i(j3) + 8007612(j|(73 - AM25Q8WV7qlql - 3iq\{92985VV7ql- 
42194i(j2'?3 - 345933v'l7^i) + 2q2{713pl - 185827Qql) + 238iVV7q3{2pl - 26Qlql) 

+ 2(?i(-1327ip§ + 6(4877^17^1 + ^eiilliqlq^ + 18207VV7q2ql - 312987iql))) 
+ Api{-24ipl{VV7qi - 2q3) - 2plp3{229qi + imq2 - 4:lVV7q3) + p3{2765VV7qt 

- 723492VV7q\ + (??(- 1466011 \/T7(/2 - 533783(^3) + 96653mqlq3 + 2909652\/T7(j|(?§ 

- 6(?f (52184 VT7(?^ - 1105425it?2(?3 - 148835VT7(j§) + 4(?i (400p§ - 550351 A/lVt/^ - 286110(/i<73 

- 873936i VT7(j2«?3 - 170A233ql) + 5A4VV7q3{pl + llSGql) + 4i(?2(620p§ + 1714059^?^)) 

+ p2(-963900i(?| + 6ip§(59\/T7(ji + 29(?3) + 42<?^(73457(?i - 9986v/T7<?3) - 18i(?i(58820(?f 
+ 36293 v/T7(?i(?3 - 121958(?i) - 7q2{50Vl7pl - 188887(?f + 24558V\7q\q3 - 2965Uqiql 
+ 4739b\/V7ql) - 9i(20825(jf + 74\9\/V7q\q3 + 95370qlql - 30634 VT7<?i(?^ + 69360^^))) 
+ 17{29097600VUql + 90(?f (20237491 \/T7(?2 + 2380085(^3) + 24(?f (13365407^17(^1 
+ 31554108i(?2(?3 + 12438441 v/T7(?i) + q\{26857pl + 12(1288345051 \/T7(?^ + 575958164^^(^3 

- 826357591 \/T7(?2i7i - 217032064(j^)) - 24(?f (-1137869911717(^2 - 714952935q\q3 
+ 1278085501 \/l7(7^^§ + 153ql{-33pl + 202589(7^) + 9iVV7q2q3{775pl + 999073ql) 

+ 29ql{273pl + 10760626(j^ ) ) + 24(7? (12944100 v'l7^2 + 1745849341(7^(73 + 25566130^ V7qlql 
+ ^(73(-1585p| + 32136222(7^) + 6iq2{-857pl + 39752817(7^)) + 16(72 (- 68397211 Vl7^f 

- 133592868(7^^3 + 8571 13651 Vl7^|^3 - 51(72(7! (13009pi + 4308412(7^) - 9lVl7^i(73(7220pi 

+ 23787879(7^) + 18(7^ (1669p§ + 41705012(7^) - 71^(p^ - 10404p|(7^ + 2255067(7^)) 

- 16(71 (8\/T7p| + p§(-306631(7^ - 52U0Vl7qlq3 + \5\mq2q3 - 4624^17(7^) 
+ 6(10068741^17(7^ - 840406051(7^(73 + 11505600VT7(7|(7§ - 2356297921(7^^^ 

- 43242781 v/T7(7l(73 + 292274371(72(73 - 2672672^17(7^))). 
For the next potential 

V7 = ^V7ql + 2401\/l4^?(72 + 330V7qiql + 935i^/Uql + 3087qlq3 + 29Aql 
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the second first integral has the following form 

h = 3030p5 + 13140p^ + 522iVUp2P3 + 96pl{120iV2p2 + llvVVps) + llp\{lQ19pl 

+ l{69576^7ql + 1775759i^\Aql + 579A7A0qlq3 + W962iVUq2q3 + 880040(?i 

+ 12ql{38327iVUq2 + 23618q3)+2AV7qi{2'i385ql - 5257ql))) + Uplp3{20iVUp3 

+ 21(37144iv/2(jf - 74257ql + 25839i\/l4<?^(?3 - 121506(?2^?3 " 58730iVUql - I2q\{3897q2 

+ 1903iVl4^3) + 6i/i(10555i V2^j^ + 6008v^(?2(?3 + 3927Qi^/lql))) + 2058p2P3(906048iVl4i/f 

- 864t?f (6577^7^2 + 6237iV2t73) + 576q\{82l3WUql + 30646(j2'?3 + 1232i Vl4tj§) 

+ 12(^1(72 (-3186975^7(71 + 58293901^2^^(73 - 5851272^7(7^(73 + 45iV2(72(79pi + 281456(7^) 
+ 4^(73 (-757p§ +403809(7^)) + (72(-7918053iVl4^f - 55828710(7|^3 - 828576iVUqlql 

- 2ql{7519pl - 35492856(7^) + i'^(>4^24ql{pl, - Mlql) - 301^(72(73 (809p§ + 1222746ql)) 
+ 12(7?(-99678iVl4(7| + 2562546(7^(73 - 2m94iVUqlql + mVUq3i41pl - 25284ql) 

+ (72(1313pi + 3566220(7^)) - 8(7?(157iV2pi + 6(695966^7(7^ - 2960791^2^^(73 

+ 604296^7(72(71 - 359758iV2(7^))) + 147pi(264p| + 2pl{2\82\.\/7q\ + 5306871^14(7^ 

+ 24i^?(6465\/T4(72 + 37031(73) + 1911084(7^^3 + 1787101 \/l4(72(7i - n(A24ql 

+ 12(71(16503^7(7^ - 392701^2(72(73 + 10598^7(7^)) + 7(-1915063849(7f + 6172(7^(107803418(71 

+ 149377659\/7(73) - 6(71(150209548(7? - 49807524a/7(7i(73 - 236319531(7^) 

+ 81^2(7^(31990100(7? + 33308910 v/7(7?(73 - 5276250(7i(7§ + 15723981^7(7^) - 12(7^(9189296(7^ 

- 3549024 V7(7?(73 - 6806016(7? + 20W792V 7qiql - 35685867(7^) + 241^2(72(898328(7? 
+ 534000^7(7? (73 - 936936(7? (?§ + 13M)220^ 7q\ql + 247m\<dqiq\ - 616665^ 7 ql) 

- 8(128528(7? - 229392 ^7(7? (73 + 720720(7? (?§ + 687456 V7(7?(7^ - 1801044(7?(7^ 

- 633276 v/7(7i (7^ - 916839(7^))) + 686 (-51 7045760^7(7? - 218437 632W\Aq\q2 

- 16128(7? ( 1231361 v/7(7^ - 7302961^2(72(73 - 116424 ^7(7^) - 9408(7? (269pi 

+ 16266621 v/T4(7^ + 8138214(7^(73 + 299376iVUq2ql + 690396(7^) - 4032(7? (526441721^14^2 

+ 375243372,71(73 - 186525361 V 14(7^ (?§ + 12iVUq2q3{227pl - 358092(7^) 

+ 924ql{5pl - 588ql) + 1813(7i(25p| + 38064^^)) + 2016(7? (p^ (271811^2(72 - 1152^7(73) 

- 56(2260588^7(71 - 25436881^2(7^(73 - 193725^/7qlql + 3441901^2^2(73 + 56448 \/7(7|)) 
+ (7i(-1803124370432iv'l4^^ - 17762238821376(7^(73 + 40916064107521 v/l4^f(7i 

+ 539101332p§(7| - 147(7|(20185577p§ - 13484416512(7^) + 181441^14^^(73 (74785p§ 
+ 62376608(7^) + 15876qlql{99587p^ + 148533504(7^) + 810lVl4^2(153p| + 374:752plql 
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- 125389824i?^)) + 36qlq2(-51113582752iVUq^ - 2482628621 76<7^ ^3 + 36273632976iVl4^|t?§ 
+ VVeUqzqli-WVpl + 112602«?^) - 392i/|(89717pi + 70945098(?^) + 126i^/Uqlq3{6M35pl 

+ 86562224(7^) +99iVl4(p| + 16268plql - 4c60992ql)) + 72qiql(-26066508160V7ql 
+ 833090912641^2(7^(73 + 18122715072\/7^|(7§ + 84^^1(73 (148769p§ - 32104800(7^) 
+ 23520iV2(72^i(-221pi + 326046(7^) + 49iV2^i(385339p§ + 67039392(7^) + 9V7(183p| 
+ 72520plql + 112482048^f )) + 96(7?(108\/7p| + 7p§ (7673471^2(7^ + 241704V7qlq3 

- 23436iV2(72(7i + 181MV7ql) - 196(65790574^7(7^ - 1283477581^2(7^(73 - \\2'^87\3\/7 q\q\ 
+ 7m624iV2qlql + 6990732^7(7^(7^ - 36991081 ^2(72(7^ - 190512^7^^))) + 2p? (-48431 V2pi 

- 9825V7pip3 + 21p2(634i^p§ + 277120iVUql + 938901^7(7^ - 18182431^2(7^(73 

+ 265062 v/7(72(73 - 6872741^2(7^ + 384,7? (2612^7(72 - 22051^2(73) + 6iqi {10499 WTiql 
+ 5946641(72(73 + 96586 v/l4(7|)) - 7p3{1040V7pl + 21(-184640(7? + 8048791^2(7^ 
+ 706404V7qlq3 - 4690981^2(72(73 + 29120V7ql + 24q\{9Q5\V2q2 + I524^7q3) 

- 96(71 (22202(7^ - 3229iVUq2q3 - 1232(7^)))) + 6pt(12319pi - 52841^14^2^3 - 7(879pi 
+ 2(12420v/7(7? - 97200iVUqlq2 - 12275301^14(7^ - 6324675(7^(73 + 6226501^14(72(73 

+ 120442(7^ + 18(71(13883^7(7^ - 470A0iV2q2q3 - 7112V7ql)))) + 6pi(4991v/2p^ 

+ \249\^7p\p3 + 7p|(-3611v^p| + 656561 ^14(7? - 421447^7(7^ + 11002951^2(7^(73 

+ 283542V7q2q3 + 6378821^2(7^ - 48(7? (1499^7(72 - 66151^2(73) + 6\qi{36241\/Uql 

+ 944721(72(73 + 12194^/Uql)) + 7p|p3(235v/7p§ + 21(50840(7? + 1543271^2(7^ 

+ 85406^7(7^^3 - 1999901 v/2(72(7i + 23828^7(7^ + (7? (-50481^2(72 - 53024^7(73) 

+ 4(7i(35881(7i + 20mVUq2q3 + 25256ql))) + 343p3 (3294720v'7(7f - 17281(7^(1117^14^2 

- 58521(73) + 2304(71(17302^7(7^ - 108431^2(72(73 - 1232^7^1) - 8(7?(1304p§ 

+ 3(310371 Vli^l + 398244(7^(73 - 7747741 Vl4(72(7i - 365344(7^)) + 12qiq2{6167733iVUqt 
+ 23478420(7^(73 + 565194:iVUqlql + 8iVUq3il39pl - 122010(7^) + 16^2(247pi + 659442(7^)) 
+ (72 ( -54123960 V7^f + 4264393141^2^^(73 + 339970176^7(7^(7^ - 573301 V2(7§(p^ + 588ql) 
+ 12\/7(72(73(6991p§ +443352(7^) - iV2ql{47317pl + 325894:296ql)) +24ql{13pl{-39iV2q2 
+ 20^7(73) + 6(832542a/7(7| - 5212551^2^^(73 + 315252^7(7^(71 - 36750iV2q2ql 
+ 25480^7(7^))) + 49p2(-1301\/2p| + 3pl{6064iVUql + 104229\/7(7f + 1193291\/2(7^(73 

- 84266v/7(72(73 - 6370iV2ql + 64(7? (1319^7(72 + 10641^/2(73) - 21,71 (10559^14(7^ 

- 714001(72(73 + 14098 v/l4(7§)) + 42 (2196481 ^2(7^ - 66545381 ^2(7^ - 12529491^7(7^(73 
+ 23141 1391 v/2(7|(7§ - 1432060^7(7^(7^ + 113195881 \/2(7^(7^ + 1078980^7(72,73 

+ 8026201^2(7^ - 6912ql{173q2 - 22iVUq3) + 2Aiq\{U8lV2ql + 299121^7(72(73 

- 18634 v/2(7§) - 4(7?(1134901(7| - 733701iVT4qlq3 + 410886(72(73 - 468861^14(7^) 

+ 121(7? (140671 v/2(7| + 9143921^7(7^(73 + 710178^ 2qlql + 3671361^7(72(73 + 177968^ 2ql) 

- 6(71 (2730321(7^ - 3671841 v/T4(7t (73 + 2040738(7i(7§ - 5460701^14(7^,73 - 660520,72(73 

+ 162681v/l4(7^)))) + 6p?(8163pt - 786\^\Aplp3 - 532lVl4p2P3 + 294p2p3(- 266921 \/2^? 
+ 49739ql - 227231 Vl4^2^3 - 9856(72^3 + 7560iVUql - 2(7? (61333(72 - 60171 Vl4(73) 
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+ qx{-(A5SeWlql + 85296v/7(j2(?3 - 2464iV2ql)) + 14p|(1132pi + 43048^7^? + 2048623iVl4(7^ 

+ 7MU74qlq3 - 86100i^/Uq2ql + 381024(7^ + 6ql(69267iVUq2 + 148666(73) 

+ 12(71 (80973^7(7^ - 2(A60iV2q2q3 - 28238^/7ql)) + 196(51p| + pl{-6100V7ql 

+ 166937iVl4(7^ + 567546^^(73 - 20055iVl4^2'73 + 59976^^ + 3^f (5361iVl4(72 + 9856^3) 

- 12(71 (8877^7(71 + 78961^2^2^73 + 910V7ql)) + 7(-258525485(7f + 3iV2(7^(40352798(7i 

+ 43391235^/7(73) - 24qj{4778174ql - 3705267V7qiq3 - 10673418(7^) + 41^2(7^(11847545(7? 
+ 7974936 V7(7?(73 - 15765939^1 (7^ + 273609V7ql) - 12(7^(1241236(7^ - 1253718^/7 qlq^ 

- 612486(7? (7§ " 78036A/7(7i(7i - 1675065qi) + 121^2(72 (284864(7^ + 126480V7q\q3 

- 481950(7f(7§ + 159012V7qlql - 572712qiqj + 13127lV7ql) - 16(54365(7^ - 25920V7qlq3 

- IWSSqfql + 2457V7qlql + ll6414qlql + 5733d^qiql - 157437ql)))). 



Vg = 27i V3990(7f + 3726Vl5(7?^2 - 456iV3996qiql - 4092Vl5ql - 1125qlq3 - 3000(7^(73 - 50ql, 
admits the following first integrals 

k = 621756pf - 171iV266pf (639p2 + 82VT5p3) - 171pl{11366pl 

+ 3381 ^p2P3 + 60(21p^ + 5(12341^3990(7? + 36^?(2159v/l5(72 + 1515(73) 

- 1^266^1 (4861 Vl5(7^ + 12780(72(73 + 4lVl5ql) - 6(3616Vl5(7| + 22665qlq3 
+ 966^15^2(73 + 35(7^)))) + 3(81377^1 + 40584Vl5plp3 

+ 75p2P3 (5586571^266(7? + 110124(7? (424(72 + 3^(73) - 2lV266(7i(1842408(7| 

+ 162352^15^2(73 - 65205(7^) - 304(72 (72658(7^ + 15753^15(72^3 + 910^^)) 

+ 750(-380805543(7f + 92341^266(7? (41202(72 - 641^(73) + 92340q\(396143ql 

- 5802^(72(73 + 235ql) - 45601^266(7? (1322406(7^ + 12177\f^q\q3 - 8505(72(7i 
+ 65^(7^) - 3420(7? (39391416(71 + 1584720^15(7^(73 + 643960qlql 

- 28980\/T5(72(7i + 3465(7|) + 304(7^(80689622(71 + 9830304^15(7^^3 

+ 9897000(7?(7§ + 650400\/l5(72(7i + 206325q^) + 241\/266(7i (234137952(7^ 

+ 18121160Vl5(7|(73 + 13755600(7^ (7§ + 406400 \/l5(7|(7i " 217350(72(73 + 875VT5ql) 

+ p§(12825i ^3990(7? + U364ql{69VT5q2 + Uqs) - 5776ql{lllVT5q2 + 196(73) 

- 21 v/266(7i (44031 v/T5(7? + 34776q2q3 - 350v/l5(7i))) + 10pi(14497p§ 

+ 15(422371^3990(7? + lUql{41736VT5q2 - 6085(73) - 21 ^266171 (233356^15(71 

+ 157860(72(73 - 4843v/T5(7|) - 76 (45750 ^15(7^ + 85660(7^(73 + 534VT5q2q3 

+ 665ql)))) + pi (631441 v/266p^ + 285961^3990^^^3 + 135p2 (3221^266^2 

+ 5(- 1374327^/15(7? + 951 v/266(7? (4706 v/l5(72 + 639(73) + 38(7i (247668 ^15(7^ 

+ 365320(72(73 - 8211v/T5(7^) - 21^266(105848^15(7? + 368340(7?(73 + 6406^15(72(73 



Potential 
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I2 = 69896714949p^ - 4410941719p^ - 12178545120 VTSp^ps 

- 5185404iv/266pf (3321p2 + llSOv/lSps) - 48735pt(9477831p^ + 6664140^15^2^3 
+ 17011399p§ + 161435514iv/3990(jf + 342^?? (31570326 v/T5(j2 + 28306975(^3) 

- 9121^266(^1 (849176 v/l5(?2 + 2324700(?2'?3 + 19775 Vl5<?§) - 4(1115132226Vl5(?| 
+ 6001374300(?i(j3 + 230353200 Vl5r?2'?3 + 208136575(?^ ) ) - 15p|(11812848011p§ 
+ 6(3637289871i V3990(?? + 57^?? (5057639934^15(^2 + 4023927275(^3) 

- 8iv/266(?i (2599271 621 v/T5(?i + 85521 17700(^2 (?3 + 1 78766525 \/T5(?3) 

- 2(54592416078v'l5(?i + 441094171900(^^(^3 + 10148787600\/l5(?2(?3 

+ 73999814225(/|))) - 420p^p3 (200059064 Vl5p§ + 75(391073751iV266^f 

+ 684(??(51068432<72 + 871817^15^3) - 2iV266<7i(1391808024i?| 

+ 204067408\/l5<72<73 - 95297139(ji) - 16(1053113326^^ + 358832133Vl5^ii/3 

+ 26186342«?2(?3 + 7756986^15^^))) - 225pi(1434130249p| 

+ 4pi (124858551871 V3990(?f - 532583214012^15(7^ - 2347905250680(7^(73 

- 28980666960Vl5(72(7i + 98949864350^^ + 171(7? (6148072746^15(72 
+ 1899245905(73) - 121^266^1 (7036440863 Vl5^| + 14574507780^2(73 

- 560965895\/l5(7§)) - 60(841855524804(7f - 584821 V266(7f (3649122(72 

+ 1662263^15^3) + 48735(71(171927564(7^ - 272077812^15(72(73 - 109530025(7§) 

- 11401 V266(7f (2587132818(7^ - 1736199048^15(7^(73 - 3090859380(72 (?! 

- 9410749^15(7^) - 1140(7? (57246559116(7^ " 28572680070^15(7^(73 

- 109359463740(7i(7i - 7143208482^15(72(71 + 1499017285(7^) 

+ 32i ^266(71 (59926048758(7^ - 24985763220^15(71(73 - 210263522175(7^(7^ 

- 20299324615^15(7^(7^ - 24297828975(72(73 + 606252325 ^15(7^) 

+ 4(766529106284,7^ - 264155752944^15(7^(73 - 9323232022800(7^(71 

- 755223062340 v/l5(7^(7l - 2978323008400(7^ (7I - 10859780400 v/l5(72(73 

- 48468636375(7^))) - 6300p2P3 (6889792 VTSpl + 15p§ (7357799971 ^266(7? 
+ 1596(7^(33673656(72 + 1541923^15(73) - 21^266,71 (2058576648(7^ 

+ 331731248^15(72(73 - 56134785(7§) - 16(1603777022(7^ + 482897263^15^2^/3 
+ 304743530(72(73 - 21530600 Vl5(7|)) - 30(54812312982^15^1 

- 877231 V266^f (305074^15(72 + 740505(73) - 1169640(7^ (1118502^15(7^ 

+ 8577605(72(73 + 98030^15^^) + 34201^266(7? (12610437^15(7^ + 705901770(7^(73 
+ 19015055V'l5^2'73 + 21915(7^) - 6840^? (273535608^15(71 - 8796044710(7^(73 

- 855723590Vl5(7|(7i + 131287940(72 (7^ + 1126965^15^^) 

+ 41^266^1 (37432776912^15(71 - 569283400800(7|(73 - 128920054050^15^1(71 

- 26461961100(7^(7^ - 407403200^15(72(73 + 1772465625(7^) + 32(24752235324^15(7^ 
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- 229256806470(jf - 110898924000 VlS^I'?! - 65789240150t?i(?i 

- 6222223575Vl5qlqi + 7013938750i?2i?3 - 242219250^15^^))) - 1125 (35195551 
+ 6p|(1651614273iV3990i?? + 171(j? (605201734 Vl5<72 + 757529455<73) 

- 16iV266(?i (466133664 Vl5(?| + 1385148660(^2 «?3 + 41393975 Vl5(j§) 

- 2(22781303074^15(7^ + 109281749540(7^(73 + 9645708800^15(72(71 

+ 879888775(7^)) + 60pi(-532302462216(7f + 175446iV266(7f (2547126(72 
+ 76525^15(73) + A38615q\{97795636ql + 6459140 Vl5^2«73 - 2853339(7§) 

- 1140i\/266(7? (5079550338(71 + 1315963350\/l5(7|(73 - 1009710576(72 (7^ 

- 7741537 v/l5(7^) - 3420(7? (28073465772(7| + 14499700170\/l5(7|(73 
+ 7760364852(7^(71 - 2099095722^15(72173 + 317762785(7|) 

+ 16iv/266(7i (180012985686(7^ + 141678458700v/T5(7|(73 + 328269338460(7^ (7^ 

- 21395932365^15(7^(7^ - 34000372350(72(73 + 289757825 ^15(7^) 

+ 4(2321951680484(7^ + 2297218340400^15(7^(73 + 10702431842160(7|(7§ 
+ 117217238580v/l5(7i(73 - 1202239869200(7^(7^ - 86811379200^15(72(73 
+ 4399443875(7^)) + 200 (1364637585276(7? (73 + 45001899iv/266(7j(166420v/l5(7| 

- 796932(72(73 + 56849v/T5(7§) + 45001899(7^(46075240^15(7^ - 215348732(7|(73 
+ 15423138Vl5^2i73 + 553855^^) - 1169641^266(7^ (6854437620 Vl5(72 

- 29454178554(7^(73 + U71374774Vl5qlql + 1180569330(72(71 - 12999235 \/l5(7|) 

- 292410(7^(134591782680^15^1 - 470046837672qjq3 - 5439697794:VT5qlql 
+ 60909452080^1 (7I + 105428700 Vl5^2'73 - 214498155^^) 

+ 684iV266(7f (5729304476700^15(71 - 14290700495904(7^(73 - 1565486188560^15(71(7^ 
+ 1995653240100(7^(7^+269861238600^15^^(7^ - 89796016800(72^^+597920425^15(7^) 
+ 2052(7? (28177936037640^15(7^ - 39924026866124(7^(73 - 11962790190192^15(7^(71 

- 3423685332600(71(7^ + 2329388397900^15(7^(71 + 1510872401600(7^(7^ 

- 55513413450^15(72(7^ + 8949071875(7^) - 288iV266(7i (5860980890620^15(7^ 

- 1807317655776(7^(73 - 2808869611056^15(7^(71 - 4526348790060(7^(7^ 
+ 328876708700^15(71(71 + 1155725294850(7^ (7^ + 10950315825 Vl5(7|(7^ 
+ 1417263750(72(73 - 206969875^15(7^) - 8(673938162327960^15(7^ 

+ 547277822406288(7^(73 - 29557137336220SVT5qlql - 943247686228260(7^ (7^ 

- 37145453338080v/T5(7^(7| + 208729566212400(7^(7^ + 12138937320150^15(7^(7^ 
+ 6956841982500(7^(7^ + 122078502000^15(72(73 + 21997219375(7^))) 

- 7980p?(-3038931iv/266p^ - 33617101^3990^^^3 + 9p2(- 19401931 V266p§ 
+ 3227221953v/l5(7? - 444031 ^266(7^(19658^15(72 + 14175(73) 

- 114(71 (148858116^15(7^ + 391208200(72(73 + 465525^15(71) 

+ 21 v/266( 181930844 v/T5(7^ + 1018181700(7^(73 + 24910050 v/T5(72(73 + 54061425(7^)) 
+ 2p3( -9905331 V3990p§ + 225(589059945^? - 10831 V266gf (136890^2 
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+ 8249^15(73) - 342(71 (8389805(7i + 1513572Vl5q2q3 + 103404^1) 

+ 2iV266(32018130^^ + 10832491^15(7^(73 + 7252524(72(71 + 9322^/^5ql)))) 

+ 855p?(216247837p| + 340148340 Vl5p^p3 + 2p|(1362509849p| 

+ 6(9274895791^3990(7? + 57(7^(1233825666^15(72 + 815729225(73) 

- 4iV266(7i(1283526583Vl5(7| + 3508986600(72(73 + 10864700 Vl5(7i) 

- 2(14147249322^15(7^ + 88856665600(7^(73 + 1646534400^15^2(73 
+ 7901396775(7^))) + 84p2P3(7399217\/T5pi + 25(98234883iv'266(7f 
+ 342(7^(21342982(72 + 860271 \/l5(73) - 2i\/266(7i (276720192(7^ 

+ 44322832^15(72(73 - 9351261(7§) - 4(834242922(7^ + 272614248^15(7^(73 

+ 96947032(72(73 - 1221759^15(7^))) + 15(53114167p| + 4p|(1073099271iv'3990(7? 

+ 171(7f (412247118v/l5(72 + 410808415(73) - 4iv/266(7i (1281131112^/15(7^ 

+ 3566728620(72(73 + 33689845 Vl5(7§) - 2(15405850098^15(7^ + 76504861120(7^(73 

+ 4369269240v/T5(72(73 - 810863525(7^)) - 20(262207731366(7^ 

- 19494iV266(7f (20214(72 +2011301^15(73) +48735(7^(153003356(7^ 

- 79345668v/T5(72(73 - 69945825(7^) - 34201^266(7^(499000974(7^ - 161887814^15(71^3 

- 328039740(72(73 - 10804817^15(7^) - 3420(7^(12226096388(71 - 2645440290^15(7^(73 

- 10784656720(7^ (7^ - 786696606^15(72(73 - 218993995(7^ )+ 48i ^266^1 (36201520638(7^ 

- 4196311620Vl5(7|(73 - 40491476550(7^ (7§ - 4439552965 v/l5(7^(7^ - 4418895600(72(73 

- 815675^15(7^) +4(1788936608636(7^ " 9012180096^15(7^(73 - 2441620201200(7|(7§ 

- 339300426660 Vl5^i(7i - 656376048600(7^(71 - 16245381600^/15(72(7^ 

- 3250847375(7^)))) - 84pi(32579496i\/266p^ + 69384500iV3990p|p3 

+ 5p|(153859518iV266pl - 110609531073^15(7? + 1711^266(7^(187871254^15(72 
+ 137031525(73) + 342(7i(1782095052Vl5(7i + 5504929400(72(73 + 68484675^15^1) 

- 21^266(5940190504^15(7^ + 45204050700(7^(73 + 887539050^15^2(73 
+ 6030275175(7^)) + 10pip3 (267620861 V3990p§ + 75( -20781023121(7? 

+ 513lV266^f (12502546(72 + 325249\/l5(73) + 342(7i (390450499(7^ + 58084012^15(72(73 

- 22450386(7§) - 21a/266(1525002138(7^ +498142221^15(71(73 +80902584(72(71 
+ 2878382\/T5(7^))) + 75p2 (90565861 v/266p| + p§(-127897609527\/l5(7? 

+ 1711^266(7? (208667026\/T5(72 + 128025855(73) + 342(7i (2125927828 \/T5(7| 

+ 5043793480(72(73 - 84854595 \/l5(7§) - 21v/266(8427582736v/T5(7i 

+ 39102611940(7^(73 + 1197339030^15(72(73 - 1363176675(7^)) 

+ 30 (- 179357568571 v/266(7^ - 555579(7^(926178(72 + 1164607^15^3) 

- 487351 ^266(7^(9819012(7^ - 6205816\/l5(72(73 - 3190425(7^) 

- 194940(7? (173644949(7^ - 63035364715(7^(73 - 107A87665q2ql - 1185437\/T5(7l) 
+ 11401 v/266(7f(2832758478(7| - 671800025^15(7^(73 - 2502926370(7^(71 

- 153606221 Vl5^2(73 + 16305030^^) + 228(7i(156415390456(7| 
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- 19923354240 Vl5^|i?3 - 180078382850(7ii?i - 18181769880 VlSt^^f?! 

- 16318512200(72 «?3 + 417338775 VlS^j^) - 41^266(135040386672(7^ 

- 752761672^15^^(73 - 198974124900(7|(7§ - 23736446770^15(7^(7^ 

- 54371797200(7^(7^ - 418389450^15^2^73 - 172155375(7^))) + 50p3 (9461481 V3990p| 
+ 5pi(-122280569379(7f + 1711^266(7^(172161162(72 + 12932861^15(73) 

+ 342(71 (1655383241(71 + 319116756^15(72(73 + 91659330(7§) - 21^266 (6445964502(7^ 

+ 2105853303^15(7^(73 + 2294558280(72(73 - 11826850\/T5(7i)) 

+ 90(-4080172176iv/3990(7f - 1666737(7^ (251098\/l5(72 + 901885(73) 

+ 48735W266qf{15997UVT5ql + 13818360(72(73 + 449075VT5(7§) 

+ 97470(7^(10497416715(7^ + 409089485(7^(73 + 17776120VT5q2ql + 7253595ql) 

+ 570iv/266(7f (110678928^15(71 - 6424828110(7^(73 - 745058085^15(7^(7^ 

- 22200210(72(73 - 5745035^15(71) + 228(7i (7552766358^15(7^ - 161449719700(7^(73 

- 40151682900\/T5(7^(7§ - 15760699275qlql + 17413200^15(72(73 + 521797500(7^) 

- 4iv/266(9011136556v/l5(7| - 114005567430(7^(73 - 56141394100VT5(7|(7i 

- 52084382850(7^ - 2311027175^15^2^3 + 3038175000^2^3 + 59134250\/l5(7|)))) 
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